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ON THERMAL CONDUCTIVITY OF COMPOSITES 



Abstract 

^Thermal conductivity of composites, made of highly 
conductive metal fibres randomly distributed in low con¬ 
ductive matrices, is investigated experimentally as well 
as theoretically* 

Experimental results are presented for a factorial 
experiment consisting of five factors at two levels.^ 
The factors are the length of the fibre, the cross 
section of the fibre, the thermal conductivity of the 
fibre, the thermal conductivity of the matrix, and the 
fibre volume ratio. 
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-/It is shown that the most meaningful factors to be 
considered for an optimum design of th^typelof compos- 
ites^in question^are, in the order of importance, the 
conductivity of the matrix, the fibre volume ratio, and 
the fibre cross section, ^ 

Earlier theoretical efforts consisted of mathemat¬ 
ical models which are highly idealized, e.g* square 
lattices, face centered cubic models (Maxwell, Rayleigh, 
Bruggeman). ^he model used in this investigation is 
that of a matrix filled with fibres, ellipsoidal in 
shape and randomly distributed, according to prescribed 
fibre volume ratios./^ It is a stochastic model,^ An 
expression of the conductivity of this composite is 
derived and results are obtained with the aid of a high¬ 
speed computer. Theoretical prediction of the stochastic 
model checks with experimental data. Also it shows that 
additional factors are involved in describing composites 
which are not covered by earlier theories. These factors 
are fibre geometry and its distribution,^ 

Moreove^r, a specific specimen formulation gives a 
thermal conductivity similar to that of stainless steels 
but a density only a quarter of that of the steel. 


Fibers (Metallic) 

Epoxy Plastics 

Manufacturing 

Methods 

Test Methods 

Statistical 

Analysis 


xi 


PART I 


INTRODUCTION 

A. Historical Review 

1. Theoretical 

Published work to date has cast light only on certain theoret¬ 
ical and experimental aspects of the thermal problems created by the 
presence of several phases of aggregates in an elastic medium. The object 
of this study is to describe experiments which serve to elucidate further 
the theory of composite thermal problems and, in particular, to show the 
significance of fiber geometry when metallic fibers are uniformly distrib¬ 
uted in low thermal conductive matrices such as epoxies. The fundamental 
equation of the thermal conductivity, and the basic relation which serves 
as a definition of that term was proposed by Fourier^ in 1822: 

Q = KAg (1) 

where Q is the amount of heat flowing per unit time through an area A, 

j 'T 

and the temperature gradient in a direction perpendicular to A is _ . 

dX 

The law states that the amount of heat flowing per second per unit area 
is proportional to the temperature gradient, and the proportionality 
constant K is the thermal conductivity. 

Equation (1) is well established for conduction in homogeneous, 
isotropic solids and is also applied for a binary mixture with two phases 
in alternate layers or laminae. In this mixture, neither phase is continuous 



is the basic equation obtained from references 2 and 3; subscript d 
denotes the discontinuous pliase; and C, the continuous phase, liquation 
(4) is found in many forms. 

In 1935, Bruggeman^ developed the model of a cubical array of 
uniform, solid spheres in a continuous medium where no limit is set on 
the amount of the phases, and no phase is necessarily completely sur¬ 
rounded by another. The equation obtained was 

K| -Kg K2 "“Kb 

Vr - ^ Vr - 

^ Kj -f2KB ^ K 2 


0 


(5) 


where the subscripts 1 and 2 refer to phase 1 and 2 respectively. 
Bruggeman also derived, from the Rayleigh-Maxwell equation, another rela¬ 
tionship applicable for dispersions of any concentrations^. Here the 
particles of conductivity are dispersed in a continuous medium of 
conductivity K2: 


(l-V, ) = 


K.-Kb 

K,-K2 



( 6 ) 


In the same work, Bruggeman also derived an equation for a mixture with 
cylindrical shapes parallel to each other and with heat current perpen¬ 
dicular to the axis of the cylinders for a wide range of concentrations: 


( Vr 


K| -Kb 
K, +Kb 


+ ( Vr >2 


Kg-Kg 
K2 -Kb 


= 0 


(7) 


Equations (1) through (7) were then generalized for any shape of the 
phases as follows^: 


( Vr 


K| Kg 
K, +FKb 


+ ( Vr )2 


K2 -Kg 
K2 + FKg 


0 


( 8 ) 


where the various values of the factor, F, result in the previously 
discussed mixture equations for the different thermal models; i.e., 

F = 0 gives equation (3) 

F = 1 gives equation ( 7 ) 

F = 2 gives equation (5) 

F = GO gives equation (2) 
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Above models are based on exact mathematical calculations, but, in 
practice, one seldom encounters such simple geometries for particles of 
varying irregular sizes. 

From 1935 and on, several intermediate formulae have been 
proposed either empirically or on the basis of some simple models. 

In 1940 and 1941, NeKrasov^ and Bogomolov*^ obtained expressions 
for heat conductivity in two extreme types of packing particles (grains): 

a. For the cubic packing: 




2 


TTKp 


0.9 - ( )m 

2.1 + (Vr)^ 


b. For the tetrahedral packing: 


( 9 ) 




3 TTKq^h 


43 4* 0. 31 P 
P - 26 


( 10 ) 


where is the thermal conductivity of the particles 

(VR)ni is the volume ratio of the matrix 

is the thermal conductivity of the air gaps 
P is the percentage of particle space (or porosity) 

In 1961, Sugawara and Yoshizawa^ developed an empirical rela¬ 
tion between experimental values of Kg and volume ratio of matrix (or 
porosity) (VR)m; i.e., 

Kb = ( I - A ) Kp-I- A Kn, 


A = 


2 " 




where 


5 


and n is an empirical exponent which they calculated as 6.5 for their 
sample. 

In 1963, Brailsford and Major^ summarized the expressions for 
the thermal conductivity of two phase media for various types of struc¬ 
tures, by assuming that each particle of phase 1 is surrounded by a 
region of material of phase 2, which is in turn surrounded by material 
having an average conductivity equal to that which has to be calculated 
(Figure 1). 

For the conductivity of a random distribution of solid spheres 
in a continuous medium, they found (by assuming that the temperature 
distribution in each region satisfies Laplace's equation = 0). 


Kg I S + Kj/Kg/^ 


+-V. 


'-X. \ 

2 + K| / Kg J 


( 12 ) 


where \f-^ is the volume fraction of phase 1. 

They also extended this model to determine the conductivity 
of a system containing three or more different phases, for example, a 
mixture in which spherical particles of material 1 or material 2 are 
embedded in material 0. The thermal conductivity of this system is given 
by: 


3K, 


3K 


KoVo+ K, V 


,-+ K, V, --r 

' ( 2K^4-K,) 2 2 (2K- -f K, ) 


K B - 


3K 


Vo+ V, 


+ V, 


(2 


3Kb I 
Ko + Kg ) J 


(13) 


( 2Ko + K,) 



Figure 1* Thermal Model Used by Brailsford and Major 
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Equation (13) assumes that phase is continuous. 

From (13) they also deduced, by letting Kq = Kg, the conductivity 
for a random two phase containing regions of both single phases in the 
correct proportions, embedded in a random mixture of the same two pliases, 
having a conductivity equal to the average value of the conductivity of the 
two phase assembly which is being calculated, i.e., 

[^(-V, - Vg )] + Kg Tsv, (2K, - Kg ) + SVg (ZKg- K , ) 




where ^2 = 1 - V^. 

In 1967, Springer and Tsai^^ introduced a thermal model which 
takes into consideration the shape of the filaments and their geometrical 
arrangement. The model is shovm in Figure 2. The basic assumptions they 
made are: 

a. Constant temperature difference between X = ± a 

b. The total heat flow per unit length along the filament, 
may be divided in three independent parts. 

c. The temperature distribution at any local point in the 
element satisfies V T = 0. 

d. The boundary conditions are: 

aty = 0, y=±b, K = 0 

dr 

at X = a, T = Ta, at x = -a, T = 

at X = 0, T = ^ ^ = constant 

2 
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e. The condition existing at the interface between the filament 
and matrix must be specified. For zero thermal resistance at the interface, 

Tf = 

.. dT . ^ 

On f ^ 

where n is the direction normal to the interface. 

With the above assumptions, the problem is then solved by 
imposing a known uniform temperature distribution (At) between x = + a 
surfaces. For symmetrical packing arrays, the following general expression 
was obtained: 


^ I - 


2a - h ) 4- (h Km /Kf ) 


Note that for a square filament (S = X = constant), a = b and 


Km " ' V 


I / Vf + B / 2 


for a cylindrical filament (S = d), a = b, hence 

-~ = I - Z^Vf/TT + -^ {tT'- ^ - 

Km B [ yi-(B^ Vf /IT 


'' I - ( B2 Vf / TT ) 
4- V B^ Vf / 77- 


where 


BS 2 


In 1968, Behrens^^ obtained thermal conductivity expressions by 

12 

using the method of analysis introduced by Max Bom into solid state 
physics, where it is commonly known as the ’’method of long waves”. It 
consists in considering plane waves, in this case thermal waves, traveling 
through the material with wavelengths which are long as compared with the 
inter-component spacings. By calculating their damping coefficients in 
the principal directions of the medium, explicit expressions for the 
average thermal conductivities can be obtained. 

For the composite with a rectangular lattice having two con¬ 
stituents (and cubic particles), in the limit, Belirens gives: 

(18) 

(P + I ) - ( P - I ) Vr 

For the composite with cubic symmetry (with spherical 
particles or inclusions), in the limit, Behrens gives: 

Kg ( P+2) + (P-l ) 2Vg ^19) 

" (P+2) - ( P-l ) Vr 

vvhere P = 

In 1968, Hashin^^ generalized the "Self-Consistent Scheme" 
or SCS method of approximation for analysis of effective properties of 
particulate composites. The basic underlying assumption is that a typical 
basic element of a heterogeneous medium, such as a single crystal in a 
polycrystal or an inclusion in a particulate composite, can be regarded as 
being embedded in an equivalent, homogeneous medium whose properties are 
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the unknowns to be calculated. This model is the same as the one introduced 
by Brailsford and Major in 1963. The general equation obtained by Hashin 

2 2 -t-c'-h/S ( I - C' ) j K^- 2 ( IH- 2C' ) + /3 ( I - 4c' ) -t- 

9(/5 - I ) V2 ] K - [ 2 ( I - C' ) + /3 ( I -+• 2C' )] = 0 

C20) 


C'= 


n *<2 

where K = — , p = - , 

K, K, 

It is noted that if C = 1, equation (20) reduces to equation (5), or 


K. - K 


K, +2Kb 


V- 


Kg Kg 
Kg + 2 Kb 


= 0 


( 21 ) 


If C = V 2 , then 


Ko = K 


V2 



( 22 ) 


This equation is also a rigorous result for a very special geometry, 
which is obtained when a volume is filled out ad infinitum with composite 
spheres of variable sizes. In each composite sphere, the volume ratio of 
the inner spherical core to matrix shell is V 2 , and the conductivities of 
core and shell are K 2 and Kj respectively. 



Hashin, in his paper, introduced C* which he believed to be a new 
geometrical parameter ’’whose value is not easily assigned,” 

3 

NOTE: If a closer look is taken at C’ = ^ , it is found tliat C’ is none 

other than the volume ratio of the particle to the matrix, and its value 
has a range of 0 to 100 percent. 

In 1789, Inger-Hauszl"^ developed one of the first techniques for 
measuring conductivities. The basic apparatus was of a comparative type; 
identical bars of different materials were coated with wax, and one end of 
the bars was heated at a constant temperature. When a steady state was 
reached, the wax was melted to different distances from the source on the 
various rods. By observing these distances, the relative thermal conduc¬ 
tivities could be determined by the equation: 



If the bars are of great length, the assumption of infinite length that is 
necessary for this equation is nearly satisfied. 

In 1861, Angstrom^^ used the impulse method of measurement of 
thermal conductivity. This method actually yields the thermal diffusivity, 
QS= so that it is necessary to have an accurate knowledge of 

specific heat and density as a function of temperature in order to determine 
the conductivity. The basic apparatus is a long bar or wire with an 
electric heater at one end. The other end extends into homogeneous surround¬ 
ings of air or an insulating medium, A steady state temperature gradient 
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is imposed on the specimen and a heat pulse varying sinusoidally vv’ith 
time is superimposed on the hotter end. This temperature disturbance is 
propagated do\vn the bar at a rate that varies with the thermal diffusivity. 

In 1867, Forbes^^ introduced the *'Forbes Bar Method". He used 
a wrought iron bar one inch in diameter and eight feet long, one end of 
which was in a bath of molten lead. The other end extended into still air. 
Heat was lost from the surface of the bar by convection and radiation. 

This introduced a longitudinal temperature gradient into the bar, and when 
a steady state was reached, tliis gradient was measured. The basic equation 
used was: 


Q 



(24) 



Following this experiment, the entire bar was heated to a higii temperature 
and allowed to cool while measurements of temperature were being made as 
a function of time. In the cooling part of the experiment, the heat 
radiated comes completely from the heat content of the bar. An accurate 
knowledge of specific heat and density is necessary to find this thermal 
energy. The energy lost per second from the bar per unit length is: 


= CA^ 


dT 
d t 


dQ 

dX 


( 26 ) 
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In 1900, Kohlrausch devised an electrical method for determining^ 
conductivity. He heated a rod by passing an electric current through 
it while holding the ends at the temperature of the surrounding. The 
method yields the ratio of thermal and electrical conductivities. The 
latter may be determined independently and accurately, therefore, the 
thermal conductivity is readily calculated. Tlie necessary measurements 
include (1) the potential drop between two points equally distant from 
the midpoint, hence at the same temperature, (2) the temperature at 
these points, and (3) the temperature, which is a maximum, at tlie middle 
of the bar. The equation used was 

K ^ ( d V / d ) _ (^27) 

^ J ( d^T/ d ) 

where 

CT is the electric conductivity (1/ohm) 

V is the electric potential (volts) 

J is the conversion factor 
18 

In 1935, Pitch used a straightforward method of measuring 
conductivity by maintaining a plate of material with one surface at a 
given temperature, then measuring the temperature of tlie opposite surface 
when a steady state of heat flow was obtained. His apparatus was unique 
in the use of a copper block for a calorimeter. Tliis apparatus is quite 
satisfactory for low conductivity materials, but for metals which have 
high heat conductivities, only a small temperature drop is encountered 
across tlie sample, leading to inaccurate conclusions. Vdien one iias a 
simple, pressed joint between two materials, as between the heat source 
and the specimen, a large and variable temperature drop is noted. This 
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drop often accounts for a large part of the temperature difference 
that is measured in connection with such a determination. This method 
is, thus, not particularly well suited for metals. 

In 1952, The American Society for Testing Materials^^ prescribed 
a method (Specs. C177-45) in which two plates of the material under test 
are placed on opposite sides of a heater. Measurements of the steady state 
temperatures of both sides of each plate gives data for the calculation 
of tlie conductivity. The assumption is made that all of the heat generated 
passes through these plates, and if a suitable guard heater is supplied 
around the edges, this condition is approximately achieved. 

In the same year (1952), Mikol^*^ modified the parallel plate 
method by using a hollow cylinder of the test material and by heating this 
cylinder at the center. The heater is in the axial hole and suitable 
guard heaters are provided at the ends. The quantity of heat flowing through 
the material is known precisely with the experimental arrangement. The 
conductivity is determined from the relationship. 


Q In (Tg / r, ) 
2 ( T| - Tg ) 


(28) 


In a commonly used experimental arrangement, the exterior of the cylinder 
is cooled by radiation either into the atmospliere or into a region of 
elevated temperature. This method is called the Radial ileat Flow Method. 

In combatting "end heat-losses," Adams and Loeb , in 1954, and 
others have used a spheroidal shaped specimen that completely surrounds 
an electric heater. This method is quite well suited for the investigation 
of the conductivity of ceramic and insulating materials. 


2 '^ 

In 1964, Goldsmid “ summarized the experimental methods and 
results of measurements on solids and fluids, as well as theories of 
heat conductions. From his discussion, it is evident that tiie modern equip¬ 
ment for measuring thermal conductivity is still of tne same type discussed 
earlier, most of the work is devoted to improve the methods of measurements 
by making them simpler, extending their range of applicability and, above 
all, by increasing their accuracy. 

The two apparati used in this study are tlie thermoconductometer 
and a comparative thermal conductivity measuring system which are described 
in detail in I’art Ill. 


B. Statement of the Problem 


No thermal data on composite materials made of high conductive 
metal fibers uniformly distributed in low conductive matrices is presently 
available in the literature* 

The property of interest in this study is the average thermal 
conductivity as it varies in relation with the following parameters: length 
of the fiber, cross section of the fiber, thermal conductivity of the fiber, 
thermal conductivity of the matrix, and fiber volume ratio* 

In addition to providing reliable thermal conductivity data, it 
is desirable to construct a theoretical model which correlates the average 
thermal conductivity and the aforementioned physical parameters. 
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PART II 


THEORY 


A. Thernal Coiiductivitv of (Composite Materials 


Tiie followinp reviev; is necessary to develop a stochastic thermal 
conductivity model. Since the model is controlled by ])robai'ilisti c lavas, 
tills study investi.i^atcs a collection of random variables from ti\e point 
of view of their interdependence and lii.-iitinp beiiavior. 

1, Probability of an Event 


X i 1 X • i '‘h 1 z., = outcomes of a random experiment 

S = tlie set of all outcomes = outcome space 
/^= a family of certain subsets S = events 
If /3 satisfies tiie following axioms: 

(C^) Tiie outcome s];ace and the erpity set are in /3 : 

s e /B j 9 e p 

If eacbi set of the finite or countable sequence , A 2 » 
Aj is in ^ , then tiicir union and intersection are in /3 : 

Ai e/3 for i = t, 2 , • • • (7 Ai e yQ 

( i) 

Ai e/S for i = I, 2 ,-• • • n Ai € /3 


(B^) If a set A is in 


tiien its complement is in 


A€ P A = S -A€ P 

then, ^ is called a Borcl field on S. Note, ayain, tliat 


S represents a certain event 


9 an impossiule event 


A an event non - 

A A) B = 9 = indicates tliat A and B are mutually exclusive (no 


elements in common) 
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kU B represents event A or B 
A/lB represents event A and B 

A set function P(A) shall be called a probability measure on ^ 
if it satisfies the following axioms: 

(Pj^) For every A in ^ , its value is a non-negative real 

number P(A)>0 for A € /5 


(P 2 ) P(S) = 1 

(ly If Aj, A 2 . ... Ay 
of mutual events, then 


is a finite or countable sequence 


PC Aj) = Z PCA; ) 
(i) ^ Ij) ^ 


For A € /3 I the designation P(A) v.'ill be called the probability of the 

event A. Note, again, that 

a. P(A) + r(A) = 1 (to prove this use P and P„) 

b. 0 < P(A) £ 1 for any event A (from a) 

c. P(Q) = 0 [from S U Q = S —P(S) + P(0) = P(S)] 

2• Probability Distribution 

l*dien a random experiment is studied, one would like to know 
the set of all its possible outcomes (outcome space), the subsets of this 
set to which probabilities are assigned (events) and, finally, the probabili¬ 
ties assigned to these su])sets. 

A probability distribution is a triple (S, yQ , P) v;here S is an 
outcome space, /3 a Borel field on S, and P a probability measure on jS , 


Iience, in order to know a probability distribution, each of its three 
components S, ^ , and P must be given. 
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3 . Continuous Probability Distributions and Probability Den sity 
Let denote the n - dimensional Euclidean space. 

A point in is an n - tuple of coordinates (X^, ... X^^) 

and can be represented in the abbreviated notation: 

(Xi, ^2> ••• ^n) ~ ^ 

In Rj^ [Line), any set defined by a < X < b, where a and b are finite 
real nur^bers, is called a closed interval. 

In R2 (plane), sets of the form, a < X^ < b, c < X2 < cl are called 
closed rectangles. In general, then, for any number of dimensions n > 1, 
sets defined by < X^ < b^, a2 < X2 < 1)2, . .. a^ < Xj^ < bj-^ will 

be called closed n- dimensional intervals. 

A probability distribution (S,y 3 , P) is called continuous when 
S is a Euclidean space R^^, yS a Borel field containing all open and close 
intervals in R^^ and P is defined as follows: 

If a function f: 

fCX^, X 2 ,...X^) = f(X) 

has the following properties: 

(Cl) f(X) > 0 for all X € Rn 

(C2) fp^ f(X) dX exists for all A£)S 

(C 3 ) /rj^ f(X] dX = 1 

Then, the probability measure P is defined as 

P(A) = J^^f(X) dX for every A 

and tlie function f satisfying Ci, C2, is called a probability density. 
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Tne integrals in C 2 and are n- dimensional integrals, and dX 
denotes the n- dimensional volume element . 

Remarks: The Borel field, j3 in Rj-^, contains all open and closed intervals, 
and all elementary geometric regions, such as interiors of closed polygons, 
curves in the plane, and solids in R 3 , with their boundaries included or 
excluded. So that 

P(A) = f f(X) dX (29) 

-'a 

ascribes a probability to any such elementary geometric region. Further¬ 
more, every bounded set A £ j3 has a volume V(A) which, for all elementary 
geometric regions, is equal to the usual volume of those regions. 

One can also show that (29) assigns probability 0 to any 
geometric region which has n- dimensional volume equal to zero. In partic¬ 
ular, no matter what the probability density f(X), probability zero will 
be assigned to single points or finite sets of points in R^, to arcs of 
curves or polygons in R 2 , to two-dimensional surfaces or pieces of sucli surfaces 
in R 3 , etc. 

NOTE: Any function, g(X) on Rn, which has the properties Cj and C 2 and 
such that 


I 



( X ) d X 9^ 0 


( 30 ) 


may be used to obtain a probability density in the form: 


f (X ) = Kg (X ) 
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for it follows from and that I > 0, hence, writing K = l/I, we have 


f(x) = ( Y ) ^ 0 ^2 


^ f(X) dX = K ^ g(X) dX exists for all A £ B 

r f(X) dX = i y g(X) dX = 1 
J 1 -. hri 


SO that f(X) satisfies C 2 , C 3 . The constant K = l/I is called the 
normalizing constant. 

4. Uniform Distribution 

Let I) be any domain in with finite positive volume; 

0 < V(D) , V < + CO 

The function 


g(X) 


1 for X € D 
0 for X € 5 


clearly has the property It also satisfies (30) since 

/ g(X) dX = r 1 dX = V 


and for any set A €/3 , we have 

f g(X) dX = r g(X) dX + / _g(X) dX = / 1 dX = V(AnD) 

AHD Ann Ano 


which is finite since it is < V(U). By using the normalizing constant 

K = 1 / V 

we obtain the probability density. 
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1 / V (D) for X € D 

0 for X 6 D 


(32) 


The continuous probability distribution on with the probability 
density (32) is known as the uniform distribution on D. 

5. Cumulative Distribution Functions 

For a one-dimensional random variable X, the function 
F(S) = P(|x:X<sj ) = P(X< S) 
is called the cumulative distribution function of X or the distribution 
function (d.f.) of X. 



b. F(S) is a non-decreasing and continuous function of S, its 
derivative exists at every point of continuity of f and, at every such 
point, we have 


d F (S ) 
dS 


f (S ) 


( 34 ) 
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c. For every pair of real numbers, < St, the jirobability 
of the event < X < S 2 1 is 


(35) 

P ( S, < X < Sg ) = F ( Sg ) - F (S, ) 


d. Also 

lim F(S) =0 
S-o -00 

lim F (S) = I 
S-^ 00 


6. Continuous n- Dimensional Random Variables 

If (X, Y) has the joint probability density f(x,y), then the 
marginal probability densities of x and y are: 

+ 00 

f (x, y ) dy 
CO 

(36) 

f ( X, y ) d X 




7* Constructing the Model 

Implementing a simulation model requires random numbers in order 
to obtain random observations from probability distributions. 

The first step is to construct the cumulative distribution 

function 

F[x) = p(x < x] 

wiiere x is the random variable involved. This can be done by writing the 
equation for this function, or grapliically plotting the function, or by 
developing a table giving the value of x for uniformly spaced values of 


I'(x) from 0 to 1. 


The second step is to generate a random decimal number between 
0 and 1, Tliis is done by obtaining a random integer number having the 
desired number of digits (including any leading zeros) and then placing a 
decimal point in front of it. 

The final step is to set p[x - x j equal to the random decimal 
number and solve for x. This value of X is the desired random observation 
from the probability function (Figure 3). 

The following procedure is recommended for locating the major 
axis of the fibers in a cylindrical shaped model (Figure 4): 


1, Determine r 


the distance from the center of the base of the 


cylinder (0 < r < R). 

2. Determine 9— the angle or location of a point on the base 

(0 < 9 < 2 TT ). 

3. Determine z —» the height of the end point within the cylinder. 

4. Introduce length of fiber i . 

5. Determine Qthe angular direction of fiber in a plane parallel 
to the xy-plane (0< CL <27r)* 

6 . DetermineCj^)-* the angular direction relative to the xy-plane 

(n< </) <2 7r). 

NOTE: 5 Random numbers are needed. 

a. To generate random observations on the base of the 
cylinder. From the uniform distribution, the joint probability density 


function 


f (x,y) = 


ttrz 


in 


elsewhere (e.w.) 


RANDOM 

OBSERVATION 


X 


Figure 3 


Illustration of Procedure for Obtaining a Random Observa 
tion From a Given Cumulative Distribution Function, 
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transforms f(x, y) into j’Cr, 0) let x = r cos 9 

y = r sin 0 


g(r, 9): f(X, y) 


^(x, y) 

a(r. 0 ) 


d{x, y) 

— 

cos 9 

sin 0 

a(r, 6) 


-r sin 9 

r cos 0 


= r 


or g(r, 9) 


IT R‘ 


'llie marginal probability densities of r and 0 are: 

zTT 

= f — ^ where 0< r < R 


and 


L(0) 


-f- 

•'o TTR 


J dr = — 
^ 27r 


w 


here 0 < 9 < ZTT 


Tlie cumulative distribution functions of f(r) and L(0) are 
r 


( r ) = / - (j r = - 

•^0 R 2 r 2 

w •‘S = 4 


where 0 < r ^ R 

where 0 < 9 < Zir 


To generate random observations, see final step in 
Constructing the Model, 

Let X = ^ , then 


r = R-yr 


or 


paragraph 7, 


G"^ (X) = r/x" 0 < X < 1 


(38) 
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Similarly, let 


-1 


y = then 0 = ITTy = F (y) 


then 

e = F"^ (u) = 27ru 0 < u < 1 (39) 

where x and u are different random numbers* 

. To generate random observations inside the cylinder * From 
the uniform distribution, the joint probability density function (for 
distributing points on a sphere) 


f(x,y,z) 


1 

J/Yu R3 


2 2 2 
in X + y + z 


0 


e.w. 


(40) 


let 


X = 


y = 


z = 


r sin 
r sin 
r cos 


in 

in 


cos 0 
sin 0 


then 

g(r, , 0) : T(x,y, z) 


^(x,y,z) 

0 ) 






sincj^ 

cos 

0 

sin (j) 

sin 

0 

cos 

4lx,. X. 

A). 


T* C 


0 

r^r\c /'T^ 

C 1 

0 

C T T1 

d Cr,0 , 



i 

-r sin 

CUo 

sin 

r cos m 

r sin (h 

S JLU 

cos 

— r Sill 

0 



0 

0 


* r^ sin 



or 


r sin cp 




3 

4 TT R. 


or sin 


g(r,(/), 0) = *{ 4 TT 


0 < Q < 2TT 
0 <<p < IT 


0 < r < R 


Tne marginal probability densities of r,ch and 9 are 


jr .zir 


3 r^ sin 


0 ^ 4 TT 


d0d<^ 


0 1 r < R 




-I I 


3 r^ sin 


0 *^0 4 77 r5 


0 sin (f) 


gdp) = 


sin 


0< (p < IT 


= f f ^ ^ ^ 

4 TT R-^ 


dr d (^ = 


1 

2 77- 
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or 


f(e) =4 


1 

ITT 

0 


in 


0 <<^< 2 


e.w. 


IT 


The cumulative distribution functions of h(r), g(c^), and f(0) are: 


IKD -J 


’’ 3 r 


dr 


R 3 


r3 


G((^) = L sin cji d. ^ ~ 


f(Q) = _i_ d 0 = _Z-_ 

Jq 2 TT 2 IT 


To generate random observations: 

for H(r) — let C = _£ 

R 



for 


for 


r = C "Mu) = R 0 ^ u ^ I (41) 

G((^) — let b = i_I_£22. — cos (p = 1 - 2b 

CL = <p = cos ^ (1 - 2b) = B ^ (u) 

a= B"^ (u) = cos"^ (1 - 2 u) 0 < u < 1 (42) 

F(0) — let a = -pff~ ® "" A ^ (0) 

(p = (u) = IWu 0 < u < 1 (43) 


where u is a random number. 


Equations (38), (39) , (41), (42) , and (43) are used with the 
length of the fiber to locate the axis of the fiber in space. 

Once this axis has been located, its position in space has 
to be established by introducing a new coordinate system. 

c. Rotation of axis in three dimensions^ From reference 23, 
x^= /jx + m^y + n^ z 

y ' = /2 ^ ^2 ^ ^2 ^ 

z' = ^ 3 X + m 3 y + n 3 z 

where ( |, m^, ni, ... n 3 ) are the cosines of the angles between Ox' 
and Ox, Ox^ and Oy, ..., Oz^ and Oz. The angles are in the direction of 
increasing to increasing x'^, where (x^, X 2 , X 3 ) corresponds to (x,y,z), 
and (x'^, x' 2 ,* ^ 3 ') corresponds to (x', y', z'). It is now necessary to 
obtain mi, n^ for i = 1,2,3 in order to locate the new coordinate axes 
X , Y , and Z '. 

Let 

= cos Qi 

= cos/Sj^ (45) 

= cosXi i ~ 1,2,3 

Recall that in paragraph 8 , Constructing the Model, ( 6 , r, z) 

were generated to obtain (Xq, Zq) which is the center of the new 

coordinate system x', y', imd z'(Figure 5a). 

(1) To find X^: Recall also that the angles oC and and the 
length of the fiber L were generated; these tliree parameters (c(j ^ L) are 
used to locate (xj, y^^, zi) , and to obtain cos Qi, cos/5i, and cos/i, 
(Figure 5b) 




Figure 5* Location of the Center and End Point of a Fiber 




Then 


= L coS(p cos Qt = Xj^ - Xq 

Ly = L coscp sin Q = yi - yo 

L 2 = L sinCj6 = Zi - zq 

^ n V 

cos Qj = —-, cospi = —— ; cos/^ = — 


(2) To find Y' (Figure 6a): The equation of the plane per¬ 
pendicular to a vector through (xq, Yq. zq) and (x^, y^, Zj) is 

(Xj - Xq) (X 2 - Xjj) + (y^ - Yq) (y2 - y^) + (z^ - z^) (z^ - z^) = 0 (48) 

is chosen so that it is always parallel to the x-y plane; hence Z 2 = 0 
and (48) becomes 


(x^ - Xq) (x 2 - X(P + (yj - Yq) (y2 - Yq) = 0 


Let the distance between (xq, Yq, Zq) and (X 2 , Yq, Z 2 ) be ^ = b; then 


(X2 - xq)^ + (y2 - Yo)^ = 


From (2) and (3): 


-(Xl - Xo^ (X2 - Xq) 


(Xi - Xq) 


for X, ^ X, 


Substituting now (51) into (50): 
(Xl - Xq)^ (X2 " Xq)^ + (72 

(Xl - Xq)-^ 




Figure 6. Rotation of Axes in Three Dimensions 
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= [(>'1 - Yq) (z2 - Zq^ ■ ^y2 - YO)] i + 

[iZi - Zq) (X2 - Xq) - (X^ - Xq) (Z2 - Zq)] J + 

[(X2 - Xo3 (Y2 - Yq) ♦ (Yi - Yq) (^2 " xq) ] k 

From preceding paragraph (2), Z 2 = Zq, hence 

z' = - (z^ - Zq) (72 - Xg) i + (zi - Zq) (X2 - Xq) j + 

[(xi - xq) (Y2 - Yq) - (Yi “ Yq) (^2 " ^q) ] k 
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where 


Z = Cj I + C2 5 + C3 ic 


Cl = X3 - xg, C2 = 73 - yg; C3 » 23 - zg 


C 0 SQ 3 = _£l—; cos^ 


3 = -cos/j = 


where 


d = 'yci + C 2 ^ + 


Paragraphs 8 c (1), (2), and (3) give the new coordinate 
system (X\ Y», Z’), with the center at (Xq, Xq, Zq) , and the following 
nine angles between (X, Y, Z) and (X*, Y’, Z*): 


cos a 2 = ~ ; COs/ 3 ^ = ~ XQ ; cosy2 = 


cos <22 = ; cos ^2 " ^ 


; COS y 7 = 0 


cos CI 3 = ; COS /33 = ^3 " ^0 ; cos y- ^3 ~ ^0 

d d d 


where 


L = [(Xl - Xg)^ + (Xl - Xg)^ + (Zj - 


i = [( 


2 2 21/2 
X 2 - Xg) + (X2 - Xo3 + (^2 " ^ 0 ) ] 


2 2 21/2 

d = [(X3 - Xg) + (X 3 - yo) + (23 - 2 o) ] 


d. Locatini 


generalize the fiber’s geometry, it is assumed that the fiber has an 


elliptical cross section. Then, from (57), L = a, > = b, and d = c are 


the distances from the center of the ellipsoid (xq, yQ, Zq) to the 
major vertex and minor vertexes respectively, 

To place an ellipsoid or fiber in space, it is necessary 
to generate 9, r, z, x, andaccording to a predetermined distribu¬ 
tion function, and then calculate the nine angles describing the 
relationship of the fiber or ellipsoid to the X, Y, Z axis. 

Tlie general equation of an ellipsoid is developed next. 

The equation of an ellipsoid centered at the origin, assuming no 
translations or rotations, is 



(58) 


where 2a is the length of the major axis, and 2b = 2c the length of 
the minor axis. 

Including rotations from equation (44), and the translations; 

x' = X - h 

y' = y - j (59) 

z' = z - k 

Tne general equation of an ellipsoid, with respect to the X, Y, Z axis, 

becomes: 


(^2 y ' ^1 °^2 ^2 

■ ■ ■ -- - ■■ - - - - + .... . + 

a2 

(60) 

(I 3 x' * y' + z')^ 

_ _ 1 


where x"^ , y\ z are given in (59), and (h, j , k) is equivalent to 
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Once the fibers are generated, they have to be placed in a 
predetermined volume (in this study, a cylindrical volume) to simulate, 
for instance, an experimental specimen. To obtain a realistic 
specimen, with uniformly distributed fibers, it is necessary to reject 
any fiber which intersects another fiber (i.e., a fiber cannot occupy 
the volume of another one) or intersects the specimen’s boundary. 

e. Method of finding the intersection of two fibers . Assuming 
two ellipsoidal fibers are given with centers at (h^, jj^, kp and (h 2 , 32 > 
k 2 ), let a > b, and a > c, then the following two conditions must be 
satisfied before the two fibers are checked for an intersection: 

[(h 2 - + Uz - 30^ * ^^2 - < 2a (61) 

and 

[(h 2 - + (j 2 - Ji)^ + ( 1^2 - > 2c (62) 

NOTE: If [(h 2 - hj)^ + (j 2 - jj)^ + (1^2 - > 2a, there is no 

possible way for the fibers to intersect, and the second fiber is kept. 

If [(h 2 - hi)^ + (j 2 - ji)^ + (k 2 - kl)^] < 2c, then the two fibers 
will intersect. Tlie second fiber is discarded, and another one is 
generated. 

When another fiber is generated, assuming that it satisfies 
conditions (61) and (62),the following procedure is used to determine if 
the fibers intersect: 

Represent the fibers as vectors and V 2 (Figure 7), where 
V]^ = ajL 1 + b]^ J + C]^ k (b3) 


and 
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^2 - ^ 2 ^ ■*' ^ 2 ^ '" 2 ^ 

Let (xj, X]^, Zj) and ?2 C ^2 » Y 2 * ^ 2 ) be the end points of 
fibers Vj and V2 respectively; and ?2> then, determine the vector 


P^P2. 


In plane geometry, nonparallel lines intersect. In solid 


geometry, however, it is possible to have nonparallel lines which do not 
intersect. In the following analysis, the perpendicular distance, d, 
between nonparallel lines is investigated. If the perpendicular distance, 
d, between the two vectors containing the fibers is zero, then and 
intersect; if d#0, but is less than the minor axis of the fibers (i.e., 
b anc c), then and V2, or the fibers, are "skew”. 

In the latter case, the coordinates of the end points of the 
perpendicular distance must be found to actually determine if the fibers 
do intersect. That is, if the coordinates of the intersections lie with¬ 
in the end points of both fibers, the fibers intersect; otherwise, there 
is no intersection. 

The perpendicular distcince between and V2 is 


d = P, P2 


V| X V2 

IV, X V2 


'Ihe coordinates of d are found next. 

Let a plane, p, pass through the point Q2 (X2, y2f ^2) > 
let the vector be normal to p; 

V4 = X V3 = + C4k (6( 


where 


( 67 ) 


Since the vector P2 Q2 perpendicular to then 

V4 . rTQ2 = 0 

hence 

Ca4i + 64) + C4k) . [(X2 - y~2 ^>"2 ' + (^2 - ^2 )^] = 0 

or 

a4 (X2 - X2') + b4 (72 - yz ') + C4 (z2 - Z2') = 0 

The equation of the plane, p, can be rewritten by multiplyiiw^ 
and transposing terms not containing the variables: 

aq X2 + b4 72 + C4 Z2 = aq X2' + bq 72' + cq Z2 

or 

aq X2 + bq 72 + cq Z2 = b (68) 

The symmetric form equations for a straight line wliich passes 
through (Xj^^ 7^, Zj^) , having the direction can be written in the 
form: 



The simultaneous solution of equations (68) and (69) yield the 
coordinates of Q^* 

The coordinates of Q2 are found as follows: 

A unit vector in the direction of (^2 

V3 _ 

— — • since Q2 has a length, d, then 

IV3I 


43 


Vs = Q2 = d + Csk ( 70 ) 

or 

q 7~Q2 = - 'V3i + (Xl' - Yl'ii * - Z2')k 

Then, the solution set (X2' , Yz' > ^2*^ obtained by 
^2' = - ^5 

72 ' = y{ - b5 ( 71 ) 

Z2' = z/ - C5 

Equations ( 64 ), (68), ( 69 ) and ( 71 ) are then used to determine 
if two fibers intersect. 

f. Method for determining whether or not a fiber is entirely 
contained within the prescribed cylindrical model . The parts to be 
considered are the top and bottom of the cylinder and the cylindrical 
envelope itself. 

( 1 ) If the z distance of the end points of tlie fiber generated 
is less than the length (L) of the model and greater than zero, the fiber 
is accepted; otherwise, it is rejected. 

( 2 ) If the X or y distance of the end points of the fiber 
generated is less than the radius (R) of the model, the fiber is accepted; 
otherwise, it is rejected. 

S. Thermal Model 

Once the fibers have been placed in the cylindrical volume, a 
method has to be devised to obtain the thermal conductivity of the 
composite model. Tlie first step is to subdivide the system into a number 
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of small but finite subvoluines, and assij^n a reference number to each, 

Tlie second step consists of writing heat-balance equations for each 
subvolume and then generalize the results. 

The simplest case of a one-dimensional heat flow, namely 
heat conduction tlirough a subvolume, will be treated. It is assumed 
that the system is exposed to a high temperature medium, i.e., a heat 
source of known and constant temperature on one side and a low-tempera¬ 
ture medium, i.e., a heat sink of known and constant temperature on tlie 
other side. Each subvolume is assumed to have an isothermal boundary 

at face A (Figure 8), (i.e., since T = T = T , heat can only flow along 

0 a 

the fictitious channel) and on insulated surface S^, (i.e., no heat is 
transferred from Si to or S ), 

If the line, z^, is dropped through the center of the subvolume, 
it will intersect the fibers at points 1 , 2 , 3 , etc.. Once these points 
are found, the channel is further simplified as shown in Figure 8, and the 
steady-state heat conduction equation can now be applied. 

Q A, AT (72) 

where ^ 


Let 


A = channel cross-sectional area 

At = - Tb 

*^Bi = Thermal conductivity of channel 1 
L = lengtl; of channel and/or model 
= length of fiber material 
= length of matrix material 

II Q 

Q = - = - r constant 

A L 




then 


(At) matrix = J} 1 L. [ (L ) i + (L ) 9 + ... + (L ) 1 


n = 1 


(At) fiber = ^ 




Q / / 

- 

Nn 


n = 1 


(L ) + ^ 

w n K 


n = 1 


(L.) 



♦ i 


which can be rewritten as: 


i— [i- 


“•»>" *r 




Let Vn be the voltune fraction of the fiber material, then 



47 



note that Vp =1 when 
1 

or 

1 

'Wi 



( 76 ) 


( 77 ) 


Knowing then ^f» ^m» can be found. 

The same procedure applies to the other channels. When , 

.—K,j are found, the overall average, Kp, of the model can be 
•^2 ^n “ 

obtained. 

a. To find fiber material in each flow channel . For various 

lines parallel to the z - axis from z = 0 to z = L, it is desired to 
determine the points of intersection between these lines and the ellip¬ 
soids. The base circle of the model, = R^, is divided into N 

parts (each part should be £ the cross sectional area of the given fiber), 
with lines going through the center of each part (Figure 8). The z's 
denote the lines. To determine the intersection of a line with an ellip¬ 
soid, we proceed as followsi assume a particular line has x and y 
coordinates of x and y. This line has to be solved with every ellipsoid 

in the direction of the heat flow (z - direction). If there is no solution, 
there is no intersection. If there is a solution 1^, z^, | Zj « Zjl is saved. 

This is done for all lines within the circle, 

b. Intersection of flow lines with fibers. Equation ( 63 ) 


is solved for z: 


where 


a = C 

b = Ex + Fy 

c = A + By^ + i)xy + Gx + Hy + I 

In this study, H = G = I = 0 

Solving now (60) for z in terms of x^ and y , the following values for 

b, and c are obtained: 

a = nf + a 2 b 2 03 ^ 


b = 2 


b^ 

{ i ^x' + nij^y') 

"3 

a"* b 

2 ( f + 1113x0] 

c = b2 

2 

c 

(t 

/ / % ^ z 

^x + m^y ) + a 

2 

2 

(t 

3x' + m3y') - a2 

a 

c 


+ r\2 a^ c2 ( ^ ) ■*■ 

(79) 

+ m2y') + 

o 2 
- c 


A computer program is provided in Appendix A to find the 
overall average thermal conductivity, when the following is given: 

a. Fiber geometry. 

b. Model geometry. 

c. Overall fiber volume ratio. 

d. Thermal conductivity of the fiber. 

e. Thermal conductivity of the matrix. 
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B. Discussion of the Planning and Analysis 

of Comparative Experiment 

1 . General Consideration in Planning the Thermal Conductivity Experiment 
of Composite Materials 

a. The nature of experimentation > An experiment has been de¬ 
fined, in the most general sense, as ”a considered course of action aimed 
at answering one or more carefully framed questions,” In the study of 
Thermal Conductivity of Composite media, however, we are concerned with 
a more restricted kind of experiment in which the experimenter does 
something to at least some of the things under study and then observes 
the effect of his action. The things under study which are being delib¬ 
erately varied in a controlled fashion may be called factors. These 
factors may be quantitative factors, such as temperature which can be 
varied along a continuous scale, or they may be qualitative factors, 
such as different compositions of matrix and fibers. The use of the prop¬ 
er experimental pattern aids in the evaluation of the factors. Many books 
have been written on the general principles of experimentation; the 
book by Wilson^"^ is used in this analysis. There are certain character¬ 
istics an experiment obviously must have in order to accomplish anything 
at all: 

( 1 ) There must be a clearly defined objective: 

(a) Choice of factors, including their range, 

(b) Clioice of experimental materials, procedure, 

and equipment. 

(c) Knowledge of what the results are applicable to. 


(2) As far as possible, effects of factors should not be 


obscured by otlier variables. The use of an appropriate experimental 
pattern helps to free the comparison of interest from the effects of 
uncontrolled variables and simplifies the analysis of the results. 

(3) As far as possible, the experiment should be free 
from bias. Some variables may be taken into account by planned group¬ 
ing; for other variables, the use of randomization is desirable. 

(4) An experiment should provide a measure of precision 
(experimental error). Replication provides the measure of precision; 
note that randomization assures validity of the measure of precision. 

(5) The experiment must have sufficient precision to 
accomplish its purpose set forth in requisite (4). Greater precision may 
be achieved by refinements of the experimental pattern technique. 

b. Experimental pattern . A common experimental pattern is the 
so-called factorial design equipment, wherein we control several factors 
and investigate their effects at each of two or more levels. If two 
levels of each factor are involved, the experimental plan consists of 
taking an observation at each of the 2^ possible combinations. 

c. Planned grouping . An important class of experimental 
patterns is characterized by planned grouping. This class is often 
called block designs. In engineering research, the tool of planned 
grouping can be used to take advantage of naturally homogeneous groupings 
in materials, machines, etc., or ’^background variables” which are not 
directly "factors” in the experiment. Statisticians have developed a 
variety of especially advantageous configurations of block designs, 
named and classified by their structure into randomized blocks, Latin 


squares, incomplete blocks, lattices, etc., with a number of sub¬ 
categories of each. 

d. Randomization . Randomization is necessary to accomplish 
requisites (3) and (4) in preceding paragraph a. In order to eliminate 
bias from the experiment, experimental variables which are not specific¬ 
ally controlled as factors, or "blocked out" by planned grouping, 
should be randomized, e.g., the allocation of specimens to treatments 
or methods should be made by some mechanical method of randomization. 

Randomization also assures valid estimates of experimental 
error and makes possible the application of statistical test of signifi¬ 
cance and the construction of confidence intervals. 

"Randomization may be thought of as insurance and, like 
insurance, may sometimes be too expensive." 

In general, we should try to think of all variables that 
could possibly affect the results, select as factors as many variables 
as can reasonably be studied, and use planned grouping where possible. 

Replication . Replication (repetition) provides the 
measure of precision, an opportunity for the effects of uncontrolled 
factors to balance out, to aid randomization as a bias-decreasing tool, 
and to spot gross errors in the measurements. 

2. Factorial and Fractional Factorial Experiments 

a. Some general remarks. In the experimental study of the 
average conductivity (K|^) of composite materials, we are interested in 
investigating the effect of length (1) , cross section (a x b) of the 
fiber, thermal conductivity of the fiber (K^), thermal conductivity of 
the matrix (K^) , volume ratio of the fiber in the matrix (Vj^) , on (K|^), 


assuming a random distribution of fibers in the matrix. 

The factors involved are then | ^ axbj K£, and 

each specific value of the factors will be called a level. 

In the past, one common experimental approach has been 
the so called "one at a time" approach. This kind of experiment would 
study the effect of varying "factors" one at a time, keeping the others 
constant. The results of such an experiment are fragmentary in the sense 
tliat we liave learned about the effect of different levels of factors at 
one factor only and vice versa. In statistical language, tliere may be an 
"interaction" effect between the two factors within the range of interest, 
and the "one at a time" procedure does not enable us to detect it. 

Factorial experiment is the name commonly applied to an 
experiment wherein we control several factors and investigate their 
effects at each of two or more levels. In the analysis of factorial experi¬ 
ments, we speak of main effects and interaction effects. Main effects of 
a given factor are always function of the average response at the various 
levels of the factor. In the case vjhere a factor has two levels, the 
main effect is the difference between tlie responses at the two levels 
averaged over all levels of the other factors. If the difference in the 
response betv/een two levels of factor A is the same regardless of the level 
of factor B (except for experimental error), we say that there is no inter¬ 
action between A and B, or that AB interaction is zero. 

If we have two levels of each of the factors A and B, then 
the AB interaction (omitting experimental error) is the difference in 
the responses of A at the second level of B minus the difference in the 
response of A at the first level of B. In general, if we have £ levels of 



53 


the factor A and ^ levels of the factor B, then the AB interaction has 

(a - 1) (b - 1) 
independent components. 

For factorial experiments with three or more factors, 
interactions can be defined similarly. For instance, the ABC interaction 
is the interaction between the factor C and the AB interaction (or 
equally between the factor B and the AC interaction, or A and the BC 
interaction). 

b. Internal estimates of errors. As in any experiment, we 
must have a measure of experimental error to use in judging the signifi¬ 
cance of the observed differences in treatments. In the larger factorial 
designs, estimates of higher-order interactions will be available. The 
usual assumption is that high-order interactions are physically impossible 
and that the estimates so labelled are actually estimates of experimental 
error. As a working rule, we often use third and higher-order interactions 
for error. The judgement of the experimenter will determine which inter¬ 
actions may reasonably be assumed to be meaningful and which may be assumed 
to be nothing more than error. These latter interactions may be combined 
to provide an internal estimate of error or a factorial experiment of 

•7 

reasonable size. For very small factorials, e.g., 2 or smaller, there 
are no estimates of high order interactions, and the experiment must be 
replicated in order to obtain an estimate of error from the experiment 
itself. 

In the case of fractional factorials, there is obviously no 
point in replication of the experiment; further experimentation would 
probably be aimed at completing the full factorial or a larger fraction 
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of the full factorial, 

identify each of the trials, in the factorial 
experiment of Kg, the following notation is adapted: 

A factor is identified by a capital letter, hence: 

= A 
axb = B 


Vr = E 

and the two levels of a factor by the subscripts: "zero^’ and ”one”. If 
we have five factors. A, B, C, D and E, then the corresponding levels of 
the factors ar Aq, Bq, ; Cq, C^; Dq, D^; Eq, E^; respectively. By 

convention, the zero subscript refers to the lower level, to the normal 
condition, or to the absence of a condition, as appropriate. A trial is 
respresented by a combination of small letters denoting the levels of the 
factors in the trial. The presence of a small letter means that the 
factor is at the level denoted by the subscript 1. The absence of a letter 
means that the factor is at the level denoted by the subscript zero. Thus 
the s>mibol ’*a*' represents the treatment combination where A is at level Aj, 
B is at level Bq, C is at level Cg, etc. The symbol "be” represents the 
treatment combination where A is at level Ag, B is at B^, C is at Cp etc. 

Conventionally, the symbol (1) represents the treatment 
combination with each factor at its zero level. In our case, Iiaving five 



combinations, and the 32 trials, are represented by 
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( 1 ) 

a, b, ab, 

c, ac, be, abc, 

d, ad, bd, abd, cd, acd, bed, abed, 

e, ae, be, abe, ee, aee, bee, abee, 

de, ade, bde, abde, ede, aede, bede, 
ab ede. 

NOTE; 1. Tlie treatment eombinations are obtained by multiplying the 
new element by all previous treatment combinations. 

2. The number of observations represents the average conductivity 
(Kg) obtained from a predetermined combination of factors. 

Figure 9 summarizes the complete factorial experiment. 
Explanation of Figure 9: 

Sample 1 is equal to Aq, Bq, Cq, Dq, Eq 
Sample c is equal to Aq, Bq, Cj, Dq, Eq 
Sample abede is equal to A^, B^, C]^, Ej^ 

The specimens used in this study are shown in Figures 10 and 11 where 

= 0.062" 

= 0.125" 

= 0.005" X 0.005" 

= 0.015" X 0.015" 

= 368 W m”^ deg"^ (Copper) 

'= 210 W m“^ deg“l (Aluminum) 

= 0.238 W m"^ deg"^ (Kpoxy) 

= 0.365 W m"^ deg“^ (Epoxy - Titanium powder) 

= 20% 
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The object is now to find the estimated main effects and interactions 
which also appear in a standard order: 

T, A, B, AB 

C, AC, BC, ABC etc... 

where 

T corresponds to the overall average effect 

A corresponds to the main effect of factor A 

AB corresponds to the Interaction of factors A and B, etc. 

d. Estimation of main effects and interactions . Yates' method 
is a systematic method for obtaining estimates of main effects and inter¬ 
action for two-level factorials. The method was originally described by 
Yates^^ and may be found in various textbooks. 

The systematic procedure for Yates' method is as follows: 

(1) Make a table with n + 2 columns. In the first column, 
list the treatment combinations in standard order. 

(2) In column 2, enter the observed response corresponding 
to each treatment combination listed in column 1. 

(3) In the top half of column 3, enter, in order, the 
sums of consecutive pairs of entries in column 2. In the bottom half of 
the column, enter, in order, the differences between the same consecutive 
pairs of entries, i.e., second entry minus first entry, fourth entry minus 
third entry, etc. 

(4) Obtain columns 4, 5..,,n+2in the same manner as 

column 3; i.e., by obtaining in each case the sums and differences of the 
pairs in the preceding column in the manner described in step 3. 
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(5) Trie entries in tlie last column (column n + 2) are 

called » corresponding to the ordered effects T, 

A, AB, etc, Hstimates of main effects and interactions are olitained 

by dividing the appropriate g by 2^ divided by is the overall 

mean. Note: Tne remaining steps of this procedure are checks on the 
computation, 

(6) ine sum of all the 2^ individual responses (colirmn 2) 
sliould ec|ual tlie total given in tne first entry of the last column 
(column n + 2) . 

(7) The sum of the squares of the individual responses 
(column 2) should equal tiie sum of the squares of the entries in the last 
column (column n + 2) divided i)y 2^. 

(8) For any main effect, ti\e entry in tlie last column 
(column n + 2) equals the sum of the resj^onses in whicii that factor is at 
its higher level minus the sum of the responses in wliicli that factor is at 
its lower level. 

e. Testing for significance of main effects and interactions , 

(Ij Choose Q , the level of significance. (See Appendix B). 

(2) If there is no available estimate of tiie variance due to 
experimental error, find tlie sum of squares of t])c g^s corresponding to 

interaction of three or raore factors. 

2 

(3) To obtain S , divide the sum of the squares obtained in 
step 2 by 2^ 1/ , v/here 1/ is the number of interaction included. In a 2^^ 
factorial, the number of third and ]iigher interactions vdll be 

Z/ = 2^ - i (n^ + n + 2) 

9 


( 80 ) 


If an independent estinate of the variation due to experimental error 
is available, use this S^. 

(4) Look up tj _ Q/2 ^ degree of freedom.If 

higher order interactions are used to obtain S^, V is the number of 
interactions included. If an independent estimate of is used, V is 
the number of degrees of freedom associated v\rith this estimate. 

(5) Compute the range V.': 


n 1/2 

= (2 ) 


1 - a/2 


(81) 


(6) For any main effect or interaction X, if the absolute 
value of gy is greater than VI, conclude that X is different from zero; 
e.g., if I g^ j > W, conclude tliat the A effect is different from 
zero. Otherwise, there is no reason to believe that X is different 
from zero. 

(7) Yates' method for obtaining estimates of main effects 
and interaction for two-level factorials can be found in Appendix C. 



C. Specimen Fabrication and Processing 


1. Introduction 

The impact and v/idespread use of fiber reinforced plastics to 
date has been largely due to development of specialized fabrication 
techniques by the manufacturer. If components with certain specified 
characteristics, such as the specimens used in this study, are desired, 
the manufacturer does not know what to do, and their suggestions are quite 
expensive. It is, then, believed that in the future composite, fabrica¬ 
tion will no longer be the sole province of the manufacturer, since it iias 
become increasingly clear that tlie designer and fabricator must work to¬ 
gether as part of an integrated effort if the full potential of fiber 
reinforced composites is to be realized for structural applications, hach 
must rely on the other for guidance. Future structural requirem.ents v;ill 
be too high to be met by the traditional empirical approach of the mold¬ 
ing shop. In turn, optimization of design is meaningless if it results in 
a structure that cannot be produced. 

2. Materials 

Each material has its o\m peculiar characteristics which deter¬ 
mine its utility in composite and which dictate how the composite must be 
processed. For instance, screw plasticization of thermosets results in a 
frozen screv;. 

The properties of each component must be understood in order 
to capitalize on its desirable properties and avoid degrading tlie result¬ 
ing composite through its misuses. The materials used in this study are: 
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a. Thermosetting resins: epoxies . Enlarged ten million times, 
the molecules of tliese resins might resemble short pieces of thread, vary¬ 
ing in length from half an inch in some of tlie liquids to several inches in 
the solids. When cured or hardened, these threads are joined togetlier at 
the ends and along the sides to form large, cross-linked structures. Each 
molecule is tied to several others, like the rope of a fish-nct or the 
filaments of a spider's web, but in an irregular, rather than neat, pattern. 
If the final products are solids, they may become soft when heated, but 
will never again liquify. The tliermosetting specimens used in this study 
are made of Epon 828 and Versamid 140 in the ratio of 60 and 40 respectively. 

"Lpon 828" is manufactured by the Shell Chemical Corporation; it 
is a liquid, aromatic epoxy in which "diepoxide 0" is the predominant con¬ 
stituent (the "0" indicates the absence of hydroxyls). 

Amine type compounds are the most widely used of all curing 
agents for epoxies, 'fliey are fast curing in room temperature, cheap, low 
in viscosity, readily miscible. Wien cured, epoxies tend to be brittle. 
Toughness, the opposite of brittleness, depends upon both tensile strengtii 
and elongation. If the stress-strain behavior of the material is measured, 
then the stress at break is the tensile strength; the strain at break is 
the elongation; and the area under the curve, the work input, is the 
toughness. Unmodified epoxies, highly cross-linked v:ith aromatic amines, 
are hard materials with high tensile strength but low elongation. At the 
otlier extreme, with moderate elongation but low tensile strength, are 
materials such as swiss cheese. Both have low areas under the stress-strain 
curve; both are brittle. Somewhere in betviieen, combining strength and 
stretch, are tough compositions needed in many industrial applications. 




Different approaches have been used to provide more stretch without 
lessenint7 the tensile strength. The ’’Versamid” resins liavc been used 
in this experiment as the curing agent for Epon 828. Tlie ’’Versamid" 
resins are viscous, broun liquids developed commercially by General 
Mills, Inc. The following advantages are found when used in epoxy formu¬ 
lations: 

(1) Viide range of compatability with epoxy resins, 

(2) Long working life. 

(3) Convenient handling without need for solvents. 

(4) No need for volatile, toxic, explosive, or 
temperature sensitive curing or accelerating agents. 

(5) Low reaction temperature, therefore, they can be cast 
in large volumes. 

(6) Low shrinkage on cure and exceptional dimensional 

stability. 

(7) Excellent resistance to meclianical impact or shock. 

(8) High flexural, compressive, and tensile strength. 

(9) Exceptional resistance to thermal shock. 

(10) Good electrical resistance, 

(11) Excellent adhesion and bonding to a wide variety of 

materials, 

(12) Good machineability. 

(13) Possibility for a wide range of properties - from hard 
and strong at one extreme, to soft and resilient on the other. 

b. Epoxy - titanium powder . This combination was chosen because 


of the excellent adhesion and bonding properties of the epoxy resin (Epon 

828 - Versamid 140) to the powder and the fibers. 
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c. Reinforcing agent s The reinforcing agent is, as the term 
implies, the primary load or heat carrying component. Consequently, for 
both efficient performance and cost effectiveness, it must be selected in 
terms of type, form, and quantity to meet the specific need. 

In tliis study, discontinuous copper and aluminum fibers were 
used to increase the thermal conductivity of the matrix resins, 

3. Manufacturing Setup and Fabrication Techniques 

Fiberglass reinforced plastics (FRP) are, by far, the best 
developed and understood of man-made fibrous composites. More types of 
composites are prepared and more resins used with glass fibers than witli 
any other reinforcing agent. Consequently, review of the fiberglass rein¬ 
forced plastics technology serves to summarize the state-of-tlie-art for 
composite fabrication and processing. 

The principal FRP fabrication methods can be conveniently 
divided into open mold and closed mold processes. Generally speaking, 
tiie open mold techniques can be characterized as yielding larger, liiglier 
performance composites, requiring more labor but less investment, having 
a slower production rate, and being highly dependent for quality on 
operator skill. 

Tlie closed mold processes, on the other hand, yield more 
reproducible composites but require much higher investment. 

In this study, a closed mold process had to be used to obtain 
specimens specified by tiie factorial experiment. Figures 12, 13, and 14 
show the ingredients and setups for manufacturing the specimens. The 
main design consideration is that of a smootii surface so that the parts 
can be easily released from the mold. It is also essential that the molds 
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Figure 13. Cutaway View of the Vacuum Casting Process 






















be rigid so that the dimensions of the part can be controlled. In this 
study, the molds were made of forged tool steel, hardened and polished, 
flie vacuum casting technique with hot-pressing v,»as used in the process of 
the specimens (Figures 12, 13, and 14). Tlie main advantage of tliis 
process is the ability to make specimens free of voids. 

The following steps should be taken to manufacture experimental 

specimens: 

a. Decide on the size of the mold and obtain the total volume 
(V^p) of the finished specimen. 

b. Weigh fibers, as specified by a predetermined volume 
ratio, per equation: 

»f ■ p, \ \ 

where 

W£ = V/eight of fibers 

= Density of fiber material 
Vj = Total volume of finished specimen 
= Fiber volume ratio 

c. Place fibers loosely in the mold. 

d. De-air matrix material by placing the proper amount of 
Rpon 828 and Versamid 140 (in different cups) in a dessicator. 

e. Mix the portions of Hpon 828 and Versamid 140 and de-air 

again. 

f. Pour this matrix on the loose fibers. 

g. Start vacuum pump to v\;et fibers. 

h. As soon as matrix shows in trap flask, stop pump, remove 
copper tubing from the mold, and set up the mold in a press. 
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i. Place heating element around the mold, start pressing 
fibers to the predetermined V^, and leave in that position for twenty 
minutes. 

j. Remove specimen from mold. 

k. If a check on the volume ratio is desired, weigh composite 
specimen and, use the follov.dng equation for V„; 



NOTE ; In this study: (1) The finished specimen v.^as cut into two samples, 
1/4** and 3/4** respectively (Figure 15 and 16), (2) The volume ratio 

of Epon 828 and Versamid 140 is 60/40 respectively. (3) The volume ratio 
of Epon 828, Versamid 140, and Titanium powder is 48-32-20 respectively. 

Hie following list shows some types of defects in the finislied 
product and the possible causes: 

Defect Possible Causes 

Blisters Cure to rapid, mold too hot, moisture in 

resin or filler. 

Pin holes, pits and voids Poor mold surface, entrapped air, in¬ 

sufficient pressure, dirty mold. 

Sink marks Cure too rapid, mold too hot, insufficient 

pressure, poor port design. 



cure 15 . Enlarged ¥iew of a 1/4 Inch Specimen 
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Possible Causes 

Resin-ricli areas caused by poor fiber 
distribution. 

Fibers content too high, poor fiber 
distribution, too much mold lubricant, 
mold too hot. 

Poor fiber distribution. 

Poor impregnation, fiber content too 
high. 




Defect 

Crazing 

Exposed Fibers 

Color variation 
Delamination 
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a. Prir-cipie^ This instrumont^ Figure 17^ was developed by 
Dr. Johan:: Schroerhn: ai Philips Zr;ntra.liaboratorium, Aachen^ Germany^ 
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will have been reached. The time taken for distilling 1 ml of liquid 
to O-mark 9 is determined v;ith a stopwatch. The thermal conductivity 
Kg is then given by 



where 

0 = Heat of vaporization for 1 ml of liquid B 

t = Time in seconds for distilling 1 ml 

T^ - Tg ^ Temperature difference in ® C which is given by the 
boiling points of the two liquids 

L “ Sample height in cm 

2 

A = Sample cross-section in cm 

R ~ Heat resistance of sample sec ® C/cal 

In this way^ a determination of the absolute thermal conduc¬ 
tivity can be made in a few minutes by a simple time measurement, When 


a number of calibrated samples are available. Kg can be determined more 
easily by a comparative measurement without knowing the exact values of 
the boiling point and the heat of vaporization of liquid B. 


Constructio n of t he apparatus , A higher boiling liquid 


A is boiled by an electric heater i at the bottom of the glass apparatus. 


The vapor passes up\^ard through pipe 3, impinges on stopper 4, flows to 
condenser 6, and returns as a liquid to the same place where it started. 

In the upper part of the apparatus 7, liquid B, at a boiling point 10-20°C 
lov/er than liquid A^ is boiled by heat from stopper 5 which has passed to 
it through sample Vapor of B passes to condenser 8 and is collected 
in graduated cylinder 9-10. 



To minimize heat loss from the upper vessel containing liquid 
B, a small boiler 2 of the same liquid B is operated in the liquid A, 
the vapor of which heats a jacket around the B liquid at the top and is 
condensed in condenser 6, and then returns to the bottom at boiler 2. 

c. Preparation and measurement . Liquid pairs can easily be 
found for the whole measuring range from 20-200®C, having suitable boil¬ 
ing points differing by 10-20®C* In Table I, some pairs are summarized 
with their boiling points and the corresponding measuring temperatures 
Tj^|, i.e., the average temperature 1/2 (T^^ + Tg). In this connection, it 
must be remembered that ether and carbon disulphide are easily flammable 
and combustible. It is important that the liquids be pure and have a 
definite boiling point. 

The heating power of the element has to be high enough so that 
liquid A boils steadily in order to get a constant heat transfer to the 
sample. To determine correct power inputs of the heating element for a 
particular liquid pair, the times are taken which are necessary to test 
a sample for different power inputs. In Figure 18, the measuring time, 
t, of the liquid pair, carbon disulphide and ether, is plotted against 
the electric power input, W. Curve II, in which the measuring time is 
independent of the power input, indicates the correct heating power. In 
curve I, the heat supply is too small, while in curve III, the boiling 
point is elevated because of the excessively high vapor pressure which 
consequently increases the temperature difference, thus decreasing the 
measuring time. In Table II, some selected values of the heating power 
are given for several liquid pairs. 

The samples, used both for calibration and measurement, are 
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TABLE I 
LIQUID PAIRS 


Liquid A 

Liquid B 

T °C 

T °C 


Methylene Chloride 

Freon 11 

41.6 

23.8 

32.7 

Carbon Disulphide 

Ether 

46,3 

34.6 

40.5 

V/ater 

Trichloroethylene 

100.0 

87.0 

93.5 

Monochlorobenzene 

Toluene 

132,0 

110.8 

121.4 

0 - D i ch 10 rob en z en e 

Monobromobenzene 

179.0 

155.6 

167.3 

SELECTED 

TABLE II 

VALUES OF THE HEATING 

POWER 




Liquid Pairs Heating Power 


Freon, 11/Methylene. Chloride 4.75 

Ether/Carbon Disulphide 4,5 

Trichloroethylene/Water 6,0 

Toluene/Monochlorobenzene 6.5 

Monobromobenzene/o-Dichlorobenzene 7,5 




• • 



cylindrical in forin and can range in diameter from 16=>18 mm and in length 
from 1--30 mm. depending on the thermal conductivity of the material« 
Samples of materials having a high thermal conductivity (pure metals) 
can be made in the form of hollov; cylinders with a wall thickness of 
about 2 miTio 

To minimize the contact resistance between sample and silver 
plates^ it is very important that the cylinder-end surfaces of the sample 
^sT efull^ ground flat and have as good a finish as can be obtained^ 
Ideallyp the surfaces of the sample should be near parallel* 

The remarkably good accuracy of measurement is due to the fact 
that the heat losses bet\/een the lower silver plate 13 and the upper 
liquid B are kept relatively low* This is achieved primarily by choosing 
the dimension of the sample so as to ensure a reasonably large heat flow 
through the sample. For instance, if the substance is a poor conductor 
of heat, the sample taken should be in the form of a thin disk. For a 
better conductor, a thicker sample or even a hollow cylinder should be 
used* ilowever, an excessively large heat flow will cause liquid B to 
boil too violently and thus decrease the accuracy of the measurement* 

The height and cross section of the cylindrical sample should therefore 
be adapted to the estimated order of its thermal conductivity. The 
measuring time should not be less than 80 seconds and not more than 1000 
secondsp while for the most exacting measurements, the time should be 
between 100 and 500 seconds (Figure 19). 

A calibration curve has to be obtained for each liquid pair 
used* Because the curve is a straight line, it is usually sufficient to 
measure only tv;o or three calibration samples* Figure 20 shows some of 




(1D3/3. 38S) M 


Nickel-Bronze 
"Remanif-Steel 
Gi Glass 
G2 Glass 
G3 Glass 

Pj^ Pressed Material 
P2 Pressed Material 
P3 Pressed Material 
P4 Pressed Material 
P<; Pressed Material 


1 (mm) 


t (see.) 


Figure 19. Calibration Diagram for an Apparatus Using the Liquid Pair, Carbon 
Disulphide/Ether as Boiling Liquids. The Known Thermal Resistance 
R of Ten Calibrated Samples is Plotted Against the Measured Time t 
Taken for 1 ml of Condensate to be Collected in Burette 14. 
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the calibration curves for several liquid pairs. They have to be 
measured anew for each apparatus. The diagrams of Figures 17 and 18 
are used here only a£ examples . 

It should be noted that for these relative measurements it 
is not necessary to know the heat of vaporization or the exact values 
of the boiling temperatures. This at once eliminates the errors that 
can be introduced into a calculated value of Kg if the boiling points 
of the liquids should differ slightly from the assumed values, owing 
to impurities, or if the ends of the samples should not attain the 
exact temperatures T^ and Tg, owing to contact resistances. In order 
to obtain reproducible results, care must be taken to ensure that these 
resistances and also the magnitude of the small, residual heat losses 
remain constant for ail measurements. Because of the protective enve¬ 
lopes, the residual heat losses are negligible in comparison with the 
great amount of heat conducted through the sample. Using the contact 
oil (graphite filled oil) and carefully ground samples, the contact 
resistances are reduced and kept constant to such an extent that their 
effect in all measurements will normally remain below the limits of 
accuracy determined by other factors. 

2. Comparative Thennal Conductivity Measuring System^ ^* 

Principle . The Comparative Thermal Conductivity Instrument, 
Figure 21, is used for determining the thermal conductivity of solid 
materials. In principle, the ability of the sample material to transport 
heat is compared to the ability of a known reference material to transport 
the same amount of heat. 


The ability to transport heat through conduction is usually de¬ 
termined by measuring the temperature difference between two points in a 













material when heat flows from one point in the direction of the other. 
The temperature difference (A T) is related to the amount of heat flow 
(Q) by Fourier's heat conduction equation; 


At= % 



(85) 


where A x is the distance between the two temperature measuring points, 
A is the area of the cross section perpendicular to the direction of 
heat flow, and K is the coefficient of thermal conductivity. 

It is quite obvious from the above equation that when the 
quantity Q/A is kept the same for two materials, one of known thermal 
conductivity (K reference) and one of unknown thermal conductivity 
sample), the respective temperature drops are related as follows; 



sample 



reference 


By measuring the temperature difference between two thermocouples a 
distance Ax apart in both reference material and test sample, the 
unknown thermal conductivity is 


y Ax j reference \ AT j sample 


( 86 ) 


The condition of equal heat flux (or heat flow density, Q/A) is obtained 
by placing a test sample and reference material of equal cross section in 
intimate contact with each other and then by holding the two between a 
heater and a heat sink. 
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In practice, a number of precautions must be taken, however, 
to ensure a constant and uniform heat flux in the sample and reference 
material. First, heat losses (or gains) along the sides will occur when¬ 
ever the temperature of the surrounding insulation does not exactly match 
those along the test stack. The fact that a temperature gradient exists 
along the test stack from the heater to the heat sink poses a difficult 
guarding problem. Second, even if the total amount of heat flow (Q) were 
kept constant, it would still be possible to obtain a non-uniform heat 
flux between the thermocouples, for example, if the disturbances caused 
by an irregular interface contact have not subsided at the location of 
the thermocouples. 

To avoid these problems, it is necessary always to have the 
surfaces of the sample and reference materials prepared to as good a sur¬ 
face finish as can be attained under the circumstances. Radial heat 
transfer is reduced by placing the test stack inside a cylindrical furnace 
with a linear temperature gradient along the length of the tube. The end 
temperatures match those at adjacent points in the stack. The space be¬ 
tween the furnace tube and the stack is filled with insulating powder. 

To obtain the best accuracy, the sample should be sandwiched 
between two identical reference materials having a thermal conductivity 
of the same order as that expected for the test material. The thermal 
conductivity is now determined as follows: 


K 


B 


2 


At j sample 



top ref¬ 
erence 


-F 



bottom 

reference 


(87) 
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1 DOUBLE-BORE ALUMINA THERMOCOUPLE 
TUBING 

2 ALUMINA BALL AND SOCKET BEADS 

3 CHROMEL-ALUMEL BARE THERMO¬ 
COUPLE WIRE 


Figure 23. Instrumentation for Pyrex and Pyroceram Reference Standards 
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The thermal conductivity of the reference materials is taken at the 
average temperature of each. In general, there are four materials which 
are used with the apparatus. These are, in the order of increasing 
thermal conductivity, Pyrex 7740, Pyroceram 9606, Inconel 702, and Armco 
iron. It is always desirable to have the temperature difference across 
the sample comparable to that across the reference materials, and thus 
the geometrical factor should be considered in relation to the choice 
of standards. Whereas the temperature gradient through a given test 
sample is determined by the output from the main heater only, the desired 
temperature level of the sample is obtained by adjusting the power into 
an auxiliary heater placed between the lower reference material and the 
heat sink. 

b. Operation and construction . After the reference materials 
and the test specimen have been instrumented according to the instructions 
given in Figures 22 and 23, the test stack can be built as shown in 
Figure 24 or Figure 25. Thermocouples are connected to the terminal 
strips on the base plate following the numbering scheme from Figure 24. 
Attention must be paid to the fact that the correct height above the heat 
sink is obtained before the lower reference standard is put in place. 

This is to ensure that thermocouples No. 1 and 2 line up properly. There¬ 
fore, an aluminum block is supplied to replace the auxiliary heater (J) 
plus insulation (K) for tests at the very lowest temperatures. A similar 
piece of aluminum is supplied to replace just the insulation (K) for 
tests at moderately low temperatures. When the furnace is lowered, the 
space surrounding the stack is filled with insulating powder, such as 
bubbled alumina. Finally, the bell jar is placed over the test section 
and the test can begin. For other details, see reference 36, 
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c. Test procedure > 

(1) Build the test stack according to instructions 
given in preceding paragraph 2a. 

(2) Lower the guard furnace. 

(3) Fasten the coolant connections in the back of 

the furnace. 

» 

(4) Fill the space inside the furnace tube with 
bubbled alumina or with some other suitable insulation. 

(5) Place the bell jar over the test section. 

(6) Check continuity of all thermocouples using the 
potentiometer or a simple ohmmeter connected across the potentiometer 
outlet on the front panel. 

(7) Turn on coolant water to the instrument and verify 
that an adequate flow passes through. (For tests using liquid nitrogen 
as the coolant, regulate the flow by hand so that only a minute amount 
of liquid nitrogen flows through the outlet tube. After a brief period 

of time a steady state will have been reached and no further flow regulation 
is necessary). 

(8) Turn all four auto transformers to zero, 

(9) Turn on the Main Power switch. 

(10) Set the red pointer of the OTC (Oven Temperature 
Control) at the desired maximum temperature and push the reset button. 

(11) Put the Temperature Reference toggle switch in the 
Temperature Reference position (this is up). 

(12) If it is desired to test in a vacuum, close the Back¬ 
fill valve and switch on the mechanical vacuum pump. 


(13) If a vacuum gauge is supplied with the instrument, 
turn on Vacuum switch. After a pressure of about 30 microns has been 
reached, the test can begin, 

(14) If it is desired to backfill the test section with 
argon or nitrogen, open the Backfill valve first and then turn off the 
vacuum pump. 

(15) If an automatic ice reference system is supplied 
with the instrument, turn on the Ice Reference switch. 

(16) Apply a small amount of power to the main, auxiliary, 
and the guard heaters. (The approximate settings of the auto transformers 
for obtaining a given test temperature level will become familiar to the 
operator after some experience has been gained with the instrument.) 

The reading (in millivolts) on the Temperature Reference meter allows a 
rapid evaluation of the main heater input, while the indicated temperature 
on the OTC may serve to provide a quick indication of the stack temperature 
level which is dependent on the auxiliary heater power. (This, of course, 
assumes that the OTC thermocouple has been placed in contact with, e.g., 
the main heater or test sample.) Power should be increased gradually at 
a rate at which the heater temperatures do not increase by more than 
about 300®F per hour. If one or two thermocouples are monitored on the 
potentiometer and the results plotted, the final equilibrium temperature 
may be roughly predicted and the power adjusted if it will be other than 
that desired. 

(17) Adjust the Lower Guard Heater to obtain TC 1 equal 


(18) Adjust the Upper Guard Heater to obtain the 
selected control thermocouple (see preceding paragraph 2b) reading equal 
to TC 7* 

(19) Wait for thermal equilibrium conditions to prevail 
throughout the test section. A sufficiently steady state has been 
reached when the temperatures do not change by more than about 1°F per 
half hour. 

(20) Record the temperatures along the test stack. Note 
i: the Ice Reference light must be cycling on and off. Note 2: the 
Temperature Reference switch must be set on TC Select before reading 

TC 4 and TC 5 on the potentiometer. 

(21) When the test is over, the vacuum may be broken by 
first opening the Backfill valve and then turning off the vacuum pump. 

If this order were reversed, the possibility exists that oil would be 
sucked from the pump into the chamber. 

d. Data reduction . In the preceding paragraph a, it was 
shown that the thermal conductivity of the sample (K^) is calculated 
from 



v/here the subscripts 

s = sample 

tr = top reference standard 

br = bottom reference standard 

This equation may be rewritten as follows: 
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The advantage is that the temperature differences may be sub¬ 
stituted by millivolt differences obtained directly from the potentiometer 
readings. Generally, this method gives more accurate results, but it 
should be kept in mind that the scheme is admissable only as long as the 
total temperature drop through the stack is sufficiently small that the 
emf versus temperature relationship of a chrome1-alumel thermocouple may 
be considered linear over this range. This is generally the case, and 
the errors in converting from emf to “F exceed those resulting from a 
possible, slight non-linearity. 

Thus, we obtain: 








(89) 


where 

Axg 




K 

s 


distance between thermocouple junctions in the sample 
in ft. 

distance between thermocouple junctions in the top 
reference in ft. 

distance between thermocouple junctions in the bottom 
reference in ft. 

sample thermal conductivity in Btu/hr-ft degF evaluated 
at a mean temperature 


TC4 ■» TC5 
2 


K^j. = thermal conductivity of the top reference in 

Btu/hr-ft degF evaluated at a mean temperature 

TC 6 + TC 7 \ 

-5- ) P 

Kjjj. = thermal conductivity of the bottom reference 

in Btu/hr-ft degF evaluated at a mean temperature 
TC 2 + TC 3 \ „ 

-3-) F 

and are obtained from the calibration curves supplied 
with the apparatus. 

Aemfg = (TC 5 - TC 4) mv 

Aemf^^ = (TC 7 - TC 6) mv 

Aemfjj^ = (TC 3 - TC 2) mv 

As a check on the accuracy of the test, one should verify that 
the total energy through the top and through the bottom reference standard 
is the same, or 
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PART IV 

DISCUSSION AND PRESENTATION OF RESULTS 

A, Materials 37, 38, 39 

The observed general characteristics of the composite specimens 
are affected by many factors, the most important of which are: 

a. Properties of the resin-hardener-additive (Ti) formulation, 

b. Properties of the fibers, 

c. Quality of the resin to fiber bond, 

d. Proportion of fiber to resin, 

e. Orientation and distribution of fibers. 

Some important aspects of performance expected from these composites are: 

a. Resistance to surface attack by water, chemicals, and 

solvents. 

b. Long-term maintenance of mechanical and thermal properties, 

c. Impermeability to liquids and gases. 

Each of the aforementioned parts could be developed in an extensive research 
program, A few important observations are presented to facilitate the 
understanding of the experimental results. Appendix B comprises most of 
the definitions of the chemical terms used in this discussion. 

As mentioned in Part II, the hardener used with Epon 828 was 
Versamid 140, This resin, made from polymerized vegetable oil acids or 
polyamid resins, is usually preferred since it has the lowest viscosity and 
highest heat distortion temperature obtainable with any Versamid. It has 
excellent adhesion to various substrates, outstanding impact resistance. 



chemical and solvent resistance, low shrinkage, and low exotherm. The 
Versamid 140 also has free primary and secondary amine groups in the 
polyamide, along with some carboxyl and a mixture of condensation products. 
IVhen mixed with epoxy resin, the amines, carboxyl, etc., are then avail¬ 
able for curing. 

Copper is one of the metals which reacts with amines, and the 
danger of both general corrosion and stress corrosion, if components are 
under stress, are present. The copper fibers used in this study were 
made from an electrolytic, tough-pitch copper sheet, the purest grade of 
copper commercially available. The minimum copper content of this material 
is 99.95 per cent, and the method of manufacture is such that no residual 
deoxidant is present (it contains only 0.04% oxygen). Because of these 
characteristics, the danger of corrosion of copper fibers, when mixed with 
the epoxy resins, is reduced. 

The aluminum fibers were made from a 99.99% commercially avail¬ 
able aluminum sheet. These fibers are highly resistant to most atmospheres 
and to a great deal of chemical agents. This resistance is due to the inert 
and protective character of the aluminum oxide film which forms on the 
metal surface. The amines are generally corrosive to aluminum, but 99.99% 

A1 is, however, resistant to amines, unless they have a very high content 
of ammonia and ammonium hydroxide. 

The Titanium powder (97.1 Ti; 0.02 Si; 0.70 N; 0.05 Fe; 0.09 H; 
0,30 Ca; 0.30 Al; 0.01 Mg; 7.6 microns is the average particle size, having 
a 99.2 - 325 mesh) is intrinsically very reactive, so that whenever the 
metal surface is exposed to air, or to any environment containing available 
oxygen, a thin, tenacious surface film of oxide is formed; this oxide film 


confers upon Titanium excellent corrosion resistance in a wide ranpe of 
corrosive media, and it behaves very favorably in contact v;ith other 
metals. 


Experimental Results 


The experimental results, obtained from the ’’Colora’’ and 
’’Comparative Apparatus”, for the factorial experiments are presented in 
Tables III through V and conveniently plotted on a semi-log scale in 
Figures 26 thru 29, The results of Yates’ factorial experiment, consist¬ 
ing of five factors at two levels, are presented in Tables VI and VII. 
Recall that the factors are the length of the fibre, the cross section 
of the fibre, the thermal conductivity of the fibre, the thermal con¬ 
ductivity of the m.atrix, and the fibre volume ratio. The analysis shows 
that the most meaningful factors to be considered for an optimum design 
of the type of composite in question are, in order of importance, the 
conductivity of the matrix, the fibre volume ratio, and the fibre cross 
section. 

In general, the results show that the aluminum fibers have a 
greater effect on the composite thermal conductivity than the copper 
fibers although the thermal conductivity of copper is almost twice that 
of aluminum. 

It would seem that there must be an added thermal resistance 
introduced between the copper fibers and the matrix. At first, it seemed 
that it could be due to a chemical reaction between the fibers and matrix, 
to an oxide film on the fibers, to the possibility that the matrix ma¬ 
terial does not adhere so well to copper as to aluminum, or a combination 
of the three. 

A chemical reaction, after reviewing the properties of the ma¬ 
terials involved, could have been possible. It is also reasonable to 
conclude that, if there was a reaction, the properties of the composite 
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Figure 27. Experimental Values of Thermal Conductivity vs. Fiber Volume 
Ratio (Al-Ep-Ti). Curve 1 (Equations 2 and 22); Curve 7 
(Equation 3). 
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Figure 29. Experimental Values of Thermal Conductivity vs. Fiber Volume 
Ratio (Cu-Ep-Ti). Curve 1 (Equations 2 and 22); Curve 7 
(Equation 3). 
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TABLE VI 


RESULTS FROM YATES' FACTORIAL EXPERIMENT ("COLORA") 


treatment 

response 

ordered 

estimated 

COMBINATION 

(KB) 

effects 

Main effect 


1 

0. 

263000E 

01 

0.994600E 

02 

0.621625E 

01 

A 

0. 

161000E 

01 

0.392000E 

01 

0.245000E 

00 

B 

0. 

104000E 

Ol 

-0.26a600E 

02 

-0.167875E 

01 

AB 

0. 

218000E 

Ol 

0.248000E 

01 

0.155000E 

00 

C 

0. 

146000E 

Ol 

0.292000E 

Ol 

0.182500E 

00 

AC 

0. 

167000E 

Ol 

0.179400E 

02 

0.112125E 

Ol 

BC 

0. 

950000E 

00 

-0.600000E 

00 

-0.375000E~01 

ABC 

0. 

191000E 

Ol 

-0.193400E 

02 

-0.i20675E 

01 

D 

0. 

453000E 

01 

0.335200E 

02 

0.209500E 

01 

AD 

0. 

298000E 

Ol 

-0.402000E 

Ol 

-0.251250E 

00 

60 

0. 

179000E 

01 

-0.180400E 

02 

-0.112750E 

01 

ABO 

0. 

229000E 

01 

-0,246000E 

01 

-0.153750E 

00 

CD 

0. 

299000E 

Ol 

-0.174000E 

Ol 

-0.108750E 

00 

ACD 

0. 

282000E 

Ol 

0.880000E 

01 

0,550000E 

00 

BCD 

0. 

292000E 

Ol 

-0.660000E 

00 

-0.412500E‘ 

-Oi 

ABCD 

0. 

241000E 

Ol 

-0.116400E 

02 

-0.727500E 

00 

E 

0. 

262000E 

Ol 

0.271000E 

02 

0.169375E 

01 

AE 

0. 

191000E 

Ol 

o.uaooooE 

Ol 

0.300000E 

00 

BE 

0. 

152000E 

01 

-0.164600E 

02 

-0.102875E 

01 

ABE 

0. 

181000E 

Ol 

-0.676000E 

01 

-0.422500E 

00 

C£ 

0. 

226000E 

01 

0.676000E 

01 

0.422500E 

00 

ACE 

0. 

453000E 

Ol 

0.151000E 

02 

0.943750E 

00 

BCE 

0. 

202000E 

Ol 

-O.aOOOOOE 

Ol 

-0.500000E 

00 

ABCE 

0. 

28500OE 

Ol 

-0.117400E 

02 

-0.733750E 

00 

DE 

0. 

100500E 

02 

0.149600E 

02 

0.935000E 

00 

ADE 

0. 

437000E 

01 

0.202000E 

01 

0.126250E 

00 

BDE 

0. 

346000E 

Ol 

-0.128000E 

02 

-0.800000E 

00 

A80E 

0. 

348000E 

Ol 

-0.600000E- 

-Ol 

-0.375000E-02 

CDE 

0. 

468000E 

Ol 

-0.378000E 

Ol 

-0.236250E 

00 

ACOE 

0. 

120500E 

02 

0.101600E 

02 

0.635000E 

00 

BCDE 

0. 

2a500oE 

Ol 

-0.506000E 

Ol 

-0.316250E 

00 

ABCDE 

0. 

282000E 

Ol 

-0.96a000E 

Ol 

-0.605000E 

00 


sample estimate 

OF VARIANCE « 2.72713200 W# 19.80450 


1 


0.994600E 

-0.268600E 

0.335200E 

0.271000E 


02 

02 

02 

02 
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TABLE VII 

RESULTS FROM YATES' FACTORIAL EXPERIMENT ("COMPARATIVE") 


treatment 

0M8INATI0N 

response 

(KB) 

ORDERED 

effects 


estimated 
Main effect 

1 

0.266000E 

01 

0.975300E 

02 

0.609562E 

01 

A 

0.164000E 

01 

0.611000E 

01 

0.381875E 

00 

B 

0.136000E 

0l 

-0.2ia700E 

02 

-0.134187E 

01 

AB 

0.262000E 

0l 

-0,997000E 

01 

-0.623125E 

00 

C 

0.141000E 

01 

0,145000E 

01 

0.906250E- 

-01 

AC 

0.213000E 

01 

0,659000E 

Ol 

0.411875E 

00 

BC 

0.1lOOOOE 

01 

-0.555000E 

Ol 

-0.346d75E 

00 

ABC 

0.210000E 

01 

-0,821000E 

Ol 

-0.513125E 

00 

D 

0.26600OE 

01 

0,266700E 

02 

0.166687E 

01 

Au 

0.259000E 

Ol 

0,569000£ 

01 

0.355625E 

00 

80 

0.200000E 

0i 

-0.941000E 

01 

-0.588125E 

00 

ABO 

0.258000E 

Ol 

-0.160300E 

02 

-0.100187E 

01 

CO 

0.261000E 

01 

-0,397000E 

Ol 

-0.248125E 

00 

ACD 

0.220000E 

01 

0.625000E 

01 

0.390625E 

00 

BCD 

0.346000E 

01 

-0,210000E 

00 

-0.131250E' 

-01 

ABCD 

0.123000E 

Ol 

. -0.555000E 

01 

-0.346875E 

00 

E 

0.239000E 

01 

0,28a300E 

02 

0.180187E 

01 

AE 

0.233000E 

Ol 

0.645000E 

Ol 

0.403125E 

00 

oE 

0.176000E 

Ol 

-0,185700E 

02 

-0.116062E 

01 

ABE 

0.160000E 

01 

-0.127500E 

02 

-0.796a75E 

00 

CE 

0.461000E 

Ol 

0,519000E 

01 

0.324375E 

00 

ACE 

0,356000E 

Ol 

0.993000E 

Ol 

0.620625E 

00 

BCE 

0.232000E 

01 

-0.661000E 

Ol 

-0.413125E 

00 

ABCE 

0.184000E 

Ol 

0,730000E 

00 

0.456250E- 

-01 

uE 

0.613000E 

Ol 

0,180500E 

02 

0.ll2ai2E 

01 

AOE 

0.760000E 

Ol 

0.138700E 

02 

0.866875E 

00 

BDE 

0.518000E 

01 

-0.915000E 

01 

-0.57ia75E 

00 

ABDE 

0.294000E 

Ol 

-0,857000E 

01 

-0.535625E 

00 

CUE 

0.326000E 

Ol 

-0.639000E 

Ol 

-0.399375E 

00 

ACDE 

0.117200E 

02 

0.155100E 

02 

0.969375E 

00 

BCDE 

0.230Q00E 

Ol 

-0.135000E 

01 

-0.843750E-01 

AoCOE 

0.314000E 

Ol 

-0,461000E 

Ol 

-0.2a8l25E 

00 


sample estimate 

OF VARIANCE a 2.6^^390160 nU 19.49995 


1 0.975300E 02 

-0.214700E 02 
0.266700E 02 
0.288300E 02 


no 


would not have changed enough to produce such high differences in thermal 
conductivities as in the case of sample ACE (Epoxy-Al fibers) which has 
a thermal conductivity almost two and a half times higher than its corre¬ 
sponding sample AE (Epoxy-Cu fibers). 

The aluminum and copper fibers usually have a thin layer of 
oxide on their surfaces. Handbook values of thermal conductivities of 
copper oxide and aluminum oxide are 1,18 and 0,68 W m“^ deg“^ respectively. 
These values, when compared with the matrices conductivities of epoxy 
and epoxy-titanium (i,e., 0,258 and 0,365 W m"^ deg“^ respectively), are 
higher; hence the possibility of an additive resistance could be disregarded. 

The bonding properties of the epoxy resin, as mentioned in Part 
II, are excellent. Some air bubbles could have been trapped during the 
specimen fabrication; however, it is believed, after inspecting the speci¬ 
mens, that their effect on the composite thermal conductivity is negligible; 
i.e., very few air holes are visible for low volume ratios and none for 
high volume ratios. 

The aluminum and copper filled specimens were made by the same 
method; hence, it was expected, at first, that the two sets of specimens 
had the same or similar fiber distribution. For the 20% fiber volume 
ratios, most of the fibers, when poured in the mold (Figure 13), filled 
the predetermined volume, evenly, and there was no need to press them. 
For the 40% fiber volume ratios, and higher, the fibers had to be pressed 
down in order to fill V^. N.B . If two specimens have identical matrices, 
identical volume ratios, identical fiber geometry, but one has aluminum 
fibers and the other copper fibers, one would expect the copper filled 
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specimens to have higher thermal conductivity than the aluminum speci¬ 
mens, since the thermal conductivity of copper is about 80^ higher than 
aluminum. 

In the 20% fiber volume ratios, the aforementioned statement 
seems to apply. Tables IV and V, Figures 26 through 29, present the 
experimental thermal conductivities and, in all cases, the 20% fiber 
volume ratios show an increase in from the aluminum specimens to the 
copper ones, in the epoxy matrix as well as the epoxy-titanium powder 
matrix. Note, for example, the samples: 


c 

-> 

Kb 

= 1.46, 

A1 

fibers 

(0.005 

X 

0.005 

X 

0.062), 

' R = 

20,8%, Ep matrix 

CD 



= 2.99, 

A1 

fibers 

(0.005 

X 

0.005 

X 

0.062), 

= 

24%, Ep-Ti matrix 

1 


’^B 

= 2.63, 

Cu 

fibers 

(0.005 

X 

0.005 

X 

0.062), 


22,4%, Ep matrix 

D 

-> 


= 4.53, 

Cu 

fibers 

(0.005 

X 

0,005 

X 

0.062), 

^R ^ 

22,3%, Ep-Ti matrix 


The increase in from C to 1 - 80% 

The increase in from CD to D = 51,5% 

The increase in Kg from C to CD = 105% 

The increase in from 1 to D = 74,5% 

From the results, it is reasonable to assume that the 20% fiber volume ratio 
specimens, having the same fiber geometry, could have similar distributions. 

In the 40% fiber volume ratios, half of the results show the 
aluminum filled specimens to be more conductive than the copper ones. This 
strange behavior could be explained if the manufacturing process and the 
mechanical properties of the fibers are considered. For this volume ratio, 
as previously indicated, the fibers had to be pressed. The copper fibers 
are stiffer than the aluminum ones; hence, when pressed, they tend to dis¬ 
tribute themselves perpendicular to the applied pressure. 
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This behavior was examined in detail by fabricating a large 
specimen, 2-1/2 inches in diameter and three inches long (Figure 30). 

From this, four specimens were machined: two, 3/4 of an inch in diameter 
and 3/4 of an inch long, to be tested in the ^^Comparative Apparatus", 
and two, 3/4 of an inch in diameter and 1/4 of an inch long, to be 
tested in the "Colora". Note that the specimens were machined in the 
direction perpendicular to the Y axis. The four specimens have a 
fiber volume ratio of 40^6 and copper fibers (0.015 X 0.015 X 0,125) in 
an epoxy matrix. As was expected, their experimental results, shora in 
Table VIII, were much higher than their corresponding ABE samples in 
Tables IV and V, This clearly indicates that specimens having copper 
fibers can be fabricated with anisotropic characteristics by just press¬ 
ing the fibers. Several private companies, universities, and governmental 
agencies, such as Wright-Patterson Air Force Base under the guidance of 
Dr. S. W. Tsai, and the University of Washington in St, Louis under the 
direction of Dr, Tolbert, are in the process of developing sophisticated 
methods (for instance, using a magnetic field) to align short fibers or 
whiskers; they feel that a successful but expensive process can be de¬ 
veloped, What would happen when the fibers or whiskers are non-magnetic? 

It is believed that the simple process found in this study could be very 
useful in some applications. In reviewing the literature on composite 
materials, no one mentions this finding. The behavior of realigning 
copper fibers under pressure can be seen if the specimens used in this 
study are closely examined. It is also noted that the realignment is 
more predominant in the ^/a ratios of 4, 8, and 25. In the epoxy matrix, 
the copper sample with the ^/a of 12,4, when compared with the corresponding 






114 


TABLE VIII 

EXPERIMENTAL RESULTS OF THE ANISOTROPIC SPECIMENS 


Fiber Vol. 


Sample 

Density 
(gm cm-3'} 

Ratio 

K 

fW m”! dee"!') 

T 

f°C) 

ABE 1 

4.16 

41 

10.5 

60.9 

ABE 3 

4.20 

42 

15.2 

59.1 

ABE 2 

4.16 

39,1 

10.23 

40.5 

ABE 4 

4.27 

40.5 

9.74 

40.5 





aluminum one, has a thermal conductivity 16% higher with the epoxy matrix 
and 115% higher with the epoxy-Ti matrix; it could be possible, then, to 
conclude that the ^/a of 12.4 gives a better or favorable fiber distribu¬ 
tion than the other ratios. 

The aluminum fibers, when pressed, showed some alignment of the 
fibers perpendicular to the applied pressure for the-^/a of 4 in the 
epoxy matrix and for the J>lz. of 8 in the epoxy-Ti matrix. (This alignment 
is not as pronounced as the one for the copper fibers.) For the other 
ratios, the aluminum fibers would bend and not relocate themselves as the 
copper ones. Note that the best fiber geometry for the aluminum fibers 
is the one with an ^/a ratio of 25 as shown in Figure 27. 

The discussion of the experimental results has been based on 
just two similar samples (the 1/4” and the 3/4”). In order to draw 
acceptable conclusions, another experiment was performed. Six 3/4 of an 
inch long and ten 1/4 of an inch long samples were chosen at random from 
a group of twenty, and a statistical analysis was carried out to test the 
repeatability or replication of results. The experimental data for this 
experiment are shown in Tables IX and X. 

With the aid of Appendix B, it is found that the mean and 
standard deviation of the conductivities, for the 3/4 of an inch samples, 
were 6.265 W m"^ deg"^ and 0.254, and for the 1/4 of an inch sample, T 
a 4.27 W m"^deg"^ and (T* 0.534. This clearly indicates that from the 
degree of consistency of the readings and reproducibility of the results, 
they are felt to be truly representative of the particular samples tested. 

Table XI shows thermal conductivity measurements on four aluminum 


specimens having 100% fibers and four fiber geometries, and four copper 
specimens having 100% fibers and four fiber geometries. 



TABLE IX 


EXPERIMENTAL 

RESULTS OF "Ai, 

Cj, Dj, E" SAMPLES 

(1/4" THICK, MEASURED 



AT 40.5'’C) 




Fiber Vol. 


Specimen 

Density 

Ratio 

K 

No. 

(gm cm‘^) 

(%) 

(W deg"^) 

1 

2.00 

39.0 

3.84 

2 

1.98 

37,5 

4.0 

3 

1.98 

37,5 

4.67 

4 

2.01 

39.0 

4.16 

5 

1.96 

35,9 

4.79 

6 

1.98 

37.5 

4.22 

7 

1.99 

38.2 

4.51 

8 

2.03 

41.4 

5.12 

9 

1.99 

38.2 

3.24 

10 

2.03 

41.4 

4.16 



117 


EXPERIMENTAL RESULTS OF "A 


TABLE X 

1» Cj, Dj, 


E" SAMPLES (3/4" THICK) 


Fiber Vol. 


Specimen 

No. 

Density 
(gm cm-^) 

Ratio 

(%) 

K 

(W deg"^) 

T 

(°C) 

1 

2.08 

35.6 

7.02 

60.8 

n 

2.07 

35.4 

5.69 

61.0 

3 

2.03 

34.0 

5.87 

60.6 

4 

2.05 

34.8 

5.81 

62.1 

5 

2.04 

34.4 

6.78 

59.8 

6 

2.06 

35.3 

6.4 

60.5 



LXPLiUMLNTAL RESULTS OF SPLCIMENS HAVING 100% FIBLR VOLUME RATIO 
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These specimens were made per equation (82), i.e., 

Vj^ is equal to 1. After the fibers were poured into the mold, a pressure 
of 29.4 tons was required to obtain the predetermined final dimensions. 
The densities of the specimens were found to be very close to the one of 
the pure metal (i.e., 2.78 gm cm"^ for aluminum and 8.92 gm cm“^ for 
copper); this indicates that the 100% fiber volume ratio was obtained. 

At this point, a very important observation is made in regard to the 
thermal conductivity measurements: the highest experimental Kg for the 
aluminum and copper specimens are 17.3 and 13.8 W m ^ deg ^ respectively; 
the discrepancy is large when these values are compared with the con¬ 
ductivities of an isotropic and homogeneous aluminum specimen (Kg = 210 
W m“^ deg"^)j or copper specimen (Kg = 368 W m"^ deg"^) as obtained from 
standard materials handbooks or by the earlier theories depicted in 
Figures 31 through 34. This experiment clearly shows the differences in 
conductivities between the pure metals and specimens made per equation 
(82), i.e., V.p Vg. In the study of such specimens, other factors 

should be considered, for example: contact resistance between fibers, 
oxide films, and other physical parameters. The purpose of this experi¬ 
ment was to clarify the 100% fiber volume ratios of the earlier theories 
and possibly open other areas of research. 


C, Tneoretical Results 

Theoretical curves obtained by other authors are shown in 
Figures 31 through 34. These curves represent the composite thermal 
conductivity versus the fiber volume ratio. Equations 9, 10, and 11, 
were not plotted since they are empirical equations and do not check 
with the results obtained in this study. Equations 16 and 17 do not 

give meaningful results; for instance, with the and K used in this 

r m 

study, a zero-division is experienced whenever volume ratios from zero 
to 100% are substituted* In obtaining the curves in Figures 31 thru 
34, it seems that it was mathematically convenient for the authors to 
assume high symmetry in representing a composite structure; for instance, 
filamentary composites were represented by square lattices, rectangular 
lattices, lamellar lattices, tetragonal lattices, cubic lattices, and 
close packed structures were represented by hexagonal lattices, face 
centered cubic etc. Equations obtained by such assumptions have rigor¬ 
ous values of composite thermal conductivity depending, in general, on 
the fiber volume ratio the thermal conductivities of the phases 

Kf), and/or some parameter determined by the mode of packing, 
or on some empirical parameter ”n”. Figures 26 through 29 definitely 
show that the factors involved in describing a composite structure are 
not only and but also the fiber geometry and its distribution. 

Curves 1 (Equations 2 and 22) and 7 (Equation 3) are plotted to 
show the upper and lov/er bounds for the thermal conductivity of composites 
as obtained by earlier theories. The stochastic thermal model developed 
in this study includes the fiber geometry and its distribution. The 
discussion of this model can be found in Part II, Part IV, and Appendix A, 
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VOLUME RATIO OF ALUMINUM FIBERS 


Figure 31. 



Thermal Conductivity vs. Fiber Volume Ratio (Al-Ep). 1. Maximum 

Conductivity (Equations 2 and 22). 2. (Equation 12). 3. 

(Equation 6). 4. Random Mixture (Equation 14). 5. (Equations 

5 and 21). 6. (Equation 7). 7. Minimum Conductivity (Equation 3). 

8. (Equation 16). 9. (Equation 18). 10. (Equation 19). 



VOLUME RATIO OF ALUMINUM FIBERS 


Figure 32. Thermal Conductivity vs. Fiber Volume Ratio (Al-Ep-Ti), 1. Max¬ 
imum Conductivity (Equations 2 and 22). 2. (Equation 12). 3. 

(Equation 6)* 4„ Random Mixture (Equation 14). 5. (Equations 

5 and 21), 6, (Equation 7). 7. Minimum Conductivity (Equation 

3). 8. (Equation 16}. 9. (Equation 18). 10, (Equation 19). 
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VOLUME RATIO OF COPPER FIBERS 


Figure 33. Thermal Conductivity vs. Fiber Volume Ratio (Cu-Ep). 1. Maximum 

Conductivity (Equations 2 and 22), 2, (Equation 12). 3. 

(Equation 6). 4. Random Mixture (Equation 14), 5. (Equations 

5 and 21). 6. (Equation 7). 7. Minimum Conductivity (Equation 

3). 8. (Equation 16), 9, (Equation 18), 10. (Equation 19). 
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Theoretical results obtained from this model are limited by 

the size of the computer, as it will be explained in section D of Part TV. 

The curves for the ^ ratios of 4, 8, 12.4 and 25 plotted in 

a 

Figures 26 through 29, have been obtained from the stochastic model. 

Note the curve ^=8: this curve is made of two parts, the section from 
a 

0 to 12.8% fiber volume ratio was obtained from the computer program in 

Appendix A and Figure 8, while from 20 to 55% was obtained from the theory 

of Appendix D and the thermal model of Figure 8. The reason for choosing 

the 8 is the following one: from the fiber volume ratio equation, i.e., 
a 

Vr = n it is noted that the number of fibers (iTi) needed to fill a 

specific volume (V-r) is n = (V^/Vf) V-g, hence for a constant VV, nwill 
increase as Vg will decrease. The Vr/V^ used in this study are listed in 
the second column of Tables XIV through XVII, Note also that the least 
amount of fibers needed to fill a predetermined volume V.j, with a prescribed 
fiber volume ratio (Vp) is the fiber with the -^/a = 8, followed by the = 

4, 25, and 12.4 respectively. 

The Mark II, G.E. computer memory bank, allowed only 500 fibers, 

or a Vr of 12.8% (see section D of Part IV) no matter what the fiber geometry, 

hence if the :^of 4 were chosen, 500 fibers would have given a fiber volume 
a 

ratio of 6.4%; the Jof 25, a V„ of 1.42%: and the :^of 12.4,a V„ of 0.71%. 

a K a 

The thermal conductivities for the aforementioned fiber volume ratios and 

geometries could have been calculated by the computer, but the answers would 

have been very close or probably the same as the matrix: for this reason 

the ii'of 8 was chosen, 
a 

Work is continuing in refining the computer program and rewriting 
the program for a larger computer such as the IBM 360. Because of the 


computer limitations, another theory was investigated. See Appendix D, 

The results from this theory are plotted in Figures 26 through 29. The 

theriaa.1 conductivities obtained from this theory follow the experimental 

behavior of the specimens, i.e., K is a function of the fiber volume 

B 

Ta-tio and fiber georaetry. The question arises: Is this theory good 

for all fiber volume ratios? By looking at the trend of the experimental 

data we can certainly make some qualitative conclusions; for example^ 

Figures 26 through 29 shov; that for tlie c::^of 8 this theory is good between 

a 

20% to 40% or perhaps higher; for the ratios of 4^ 12.4, and 25 the 

a 

theory could be useful for fiber volume ratios higher than 40%. The 
computer program of Appendix A is believed to be correct; the limited 
ds'ta obtained from it seems to be directed toward the experimental values, 
only a larger computer could ascertain its validity. 



D, Computer Program 


The computer program is written for the Mark II G.E. time 
sharing system. The program is sufficiently versatile so that it can 
be run by anyone. To gain access to the computer, one depresses the 
ORIG button on the teletypewriter console and dials the telephone 
number of the central computer. The computer then responds with a 
series of statements. The underlined portions are user supplied while 
the remainder are computer responses. 


USER NUMBER - - - - SBT45111. TUMRDE 

PROJECT ID - - - - FIBERS 

SYSTEM .TOR 

NEW OR OLD.0^ 

OLD FILE NAME - - - COMP 1 
READY 

RUN 

COMP 1 09.38 03/28/69 

FORT-C-F DATED 10 MAR 1969 
SAMPLE SIZE INCREMENT? 10 


(The number 10 is typed here to generate 10 fibers into the specified model 
volume). 

FOR SAMPLE SIZE = 10 DISCARDED FIBERS 0 

(This statement tells tlie user that no fibers intersected any other fiber 
when introduced in the matrix. NOTE: The fibers intersecting the model's 
surface are rejected, but are not counted in "DISCARDED FIBERS".) 

MORE FIBERS? 0_ 

(The user types 0 or 1 if his answer is "yes" or "no" respectively.) 





SAMPLE SIZE INCREMENT? 20 


FOR SAJ-1PLE SIZE = 20 DISCARDED FIBERS 2 

(Here 2 new fibers intersected with one of the previous fibers; there¬ 
fore they were discarded and tv?o ne'w ones were generated in their place,) 
MORE FIBERS? 

FIBER DENSITY CHECK? 0 

(Again, the user types 0 or 1 for yes or no. If the answer is yes, the 
user will now perform the heat conduction problem. At first the model 
is subdivided in a number of heat channels, hence: 

VALUES FOR DY AND DX? .015, .015 

(Here the user states the size of the integrating unit or heat channel 
area that is desired to be used. In this case, an area of .000225 is 
used to estimate the fiber density, since the fiber used in the analysis 
has a cross section of 0.015 X 0,015 square inches.) 

FIBER HITS = 1.0300E + 02 

MATRIX HITS = 1.8730E + 03 

SUMKB = 2.9037E - 01 
SUMKM = 4,4578E + 00 
AVERAGE KB = 2.4029E - 03 

(FIBER HITS represent the number of heat channels where fibrous material 
wfas located. I^Al'RIX HITS represent the number of heat channels where 
fibrous material was not located. SUMKB represents the sum of the thermal 
conductivities of the FIBER HITS. SUMKM represents the sum of the thermal 
conductivities of the MATRIX HITS. Note that (MATRIX HITS) times (K^^) 
should be equal to SUMKM, in this case K^ = 0.0023S V.' cm'^ deg’^. AVERAGE 
Kp is the total thermal conductivity of the matrix with the specified 


number of fibers or the specified fiber volume ratio,) 

FIBER DENSITY CHECK? 

(If the answer were 0, the computer would again respond with a request 
for: "VALUES FOR DY AND DX?", thus enabling the user to change the in¬ 
tegrating interval if desired.) 

MORE FIBERS? STOP 

(The program is terminated with this response. If more fibers were 
desired and the proper response typed, the user could go through the 
aforementioned cycle again. NOTE: The program could be terminated at 
any response by typing STOP or typing ANY NUMBER GREATER THAN ONE AT 
"MORE FIBERS?"). 

This computer program was unable to generate more than 500 
fibers since the program required a considerable amount of computer 
memory, and the Mark II G.E, Time sharing system was unable to furnish 
the proper amount. 

Tables XII and XIII present computer results for the average 
thermal conductivity of a model having a geometry similar to the "ABC" 
specimen used in the "COLORA" instruments. These results were plotted 
in Figure 26 for comparison with the experimental ones. Some of the 
disadvantages of this program are the limitation on the fiber volume 
ratio and the cost of the computer time; for instance, the running time 
for Table XIII was 1664 units, the terminal time was 146 units. One 
unit costs 30 cents; hence 1810 units cost $543.00! (Note that one 
unit is equivalent to 4/3 of a second), or Table XII cost $2172.00. 
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TABLE XII 


No. of 
Fibers 

COMPUTER RESULTS FOR 

FIBER: 0.015" 

X 0.015" X 

0.125” 

(.h 


Kg (W deg"^) 


Al-Ep 

Al-Ep-Ti 

Cu-Ep 

Cu-Ep-Ti 

0 

0 

0.238 

0.365 

0.238 

0.365 

20 

0.51 

0.240 

0.367 

0.240 

0.367 

100 

2.56 

0.252 

0.370 

0.252 

0.370 

200 

5.12 

0.376 

0.441 

0.456 

0.520 

300 

7.68 

0.503 

0.549 

0.663 

0.709 

350 

8.98 

0.517 

0.557 

0.677 

0.714 

500 

12.80 

0.923 

0.945 

1.334 

1.356 


TABLE XIII 
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EXAMPLE OF COMPUTER PROGRAM OUTPUT 

COMPl 15:17 04/11/69 

FORT-C-F DATED 10 MAR 69 

SAMPLE SIZE INCREMENT7100 

FOR SAMPLE SIZE = 100 DISCARDED FIBERS 14 

MORE FIBERS?! 

FIBER DENSITY CHECK70 

FIBER CONDUCTIVITY CONSTANT72.1 

VALUES FOR DY AND DX?.015,.015 

FIBER HITS = 5.7700E+02 

MATRIX HITS = 1.3990E+03 

SUMKB = 1.6649E+00 SUMKM = 3.3296E+00 

AVERAGE KB = 2.5276E-03 

FIBER DENSITY CHECK?1 

MORE FIBERS70 

SAMPLE SIZE INCREMENT7400 

FOR SAMPLE SIZE = 500 DISCARDED FIBERS 843 

MORE FIBERS?! 

FIBER DENSITY CHECK70 

FIBER CONDUCTIVITY CONSTANT72.1 

VALUES FOR DY AND DX?.015,.015 

FIBER HITS = 1.6250E+03 

MATRIX HITS = 3.5100E-I-02 

SUMKB = 1.7414E+01 SUMKM = 8.3538E-01 

AVERAGE KB = 9,2354E-03 

FIBER DENSITY CHECK70 

FIBER CONDUCTIVITY CONSTANT73.68 

VALUES FOR DY AND DX?.015,.015 

FIBER HITS = 1.6250E+03 

MATRIX HITS = 3.5100E+02 

SUMKB = 2.5522E+01 SUMKM = 8.3538E-.01 

AVERAGE KB = 1.3339E-02 

FIBER DENSITY CHECK?1 

MORE FIBERS710 

PROGRAM STOP AT 1520 
USED 1373.28 UNITS 



E, Distribution of Fibei- s in a 5-Diir.e nsj^onal S pace 


cofTiir^uter m 
Appendix D 
betu'een an 
ecuivalent 


The Mark II computer was unable to furnish the proper 

S'-ory bank for fiber volume ratio higher than 12,8, In 

C 

an analysis is presented to detennirie the mean distance^ O ^ 
individual fiber and its nearest neighbor* This distance is 
to a. mean amount of matrix material betv/een a fiber and its 


nearest neighbor in a collection of fibers randomly distributed in a 3- 
dimensional space* It is noted that in the expression for S , the fiber 
vciuma ratio and the fiber geometry play a very important role in the 
later determination of the composite thermal conductivity* Results 
obtained from this approach^ for the 20 and 40 percent fiber volume ratiOj, 
seem to have the same behavior as the experimental results and, in some 
cases,, for the aluminum fibers with an J/a of 8, the match is very good. 
Tables XIV througii XVII show theoretical results obtained from this theory 
for volume ratios of 20, 30, 40, 50 and 60% for the different geometries 
and materials used in this study, NOTE: This model does not work for 
volume ratios of 0 and 100% (see the expression for (S in Appendix D), 
Probably more reliable theoretical models could be developed to piredict 
the thermal conductivity of the type of composites used in this study. 

The experimental work indicates the general behavior of the thermal con¬ 
ductivity of composites as a function of Vp^, the fiber geometry, 

and its distribution; this data should be very helpful in developing 
better and more reliable models* Because of the many variables encountered 
it seems almost impossible to approach this problem from a deterministic 
point of viev;. In this study an attempt v/as made to develop a statistical 
cr stochastic model, the theoretical and experimental results seem to indi¬ 
cate the feasibility of the statistical and probabilistic approach. 
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TABLE XIV 

THERMAL CONDUCTIVITY OF FIBERS RANDOMLY DISTRIBUTED IN A 3-D SPACE 
FIBER: 0.005" X 0,005" X 0.062" 


Fiber Vol. 
Ratio 
(%) 

Vr 

Vf 



*Kg (W 

deg"^) 


(cm“3) 

L 

Al-Ep 

Al-Ep-Ti 

Cu-Ep 

Cu-Ep-Ti 

20 

7874 

0.043 

5.29 

8.02 

5.35 

8.15 

30 

11811 

0.038 

6.04 

9.13 

6.11 

9.30 

40 

15748 

0.034 

6.63 

10.00 

6.72 

10.21 

50 

19685 

0.032 

7.12 

10.74 

7.23 

10.97 

60 

23622 

0.030 

7.55 

11.37 

7.67 

11.64 


*Kj, = Thermal Conductivity 


TABLE XV 

TIlERf-IAL CONDUCTIVITY OF FIBERS RANDOMLY DISTRIBUTED IN A 3-D SPACE 
FIBER: 0.005" X 0.005" X 0.125" _ 


Fiber Vol. § *Kg (W m"^ deg"^) 

Ratio Vf “T __ 



(cni“3) 

JU 

Al-Ep 

Al-Ep-Ti 

Cu-Ep 

Cu-Ep-Ti 

20 

3906 

0.055 

4.21 

6.40 

4.25 

6.48 

30 

5859 

0.048 

4.81 

7.29 

4.86 

7.40 

40 

7812 

0.044 

5,28 

8.00 

5.34 

8.13 

50 

9765 

0.041 

5.68 

8.59 

5.74 

8.74 

60 

11719 

0.038 

6.03 

9,11 

6.10 

9.27 


*Kg = Thermal Conductivity 
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TABLE XVI 


THERMAL CONDUCTIVITY OF FIBERS R.ANDOMLY DISTRIBUTED IN A 3-D SPACE 




FIBER; 

0,015” X 0,015” X 0,062” 

Fiber Vol. 

Vr 

s 

(V\ dep'^l) 

Ratio 

Vf- 


(-0 

L 

Al-Ep Al-F,p-Ti Cu-Ep Cu-Ep-Ti 




20 

874e5 

0.092 

2.58 

3.93 

2.59 

3,96 

30 

1311,8 

0.079 

2,95 

4.49 

2.96 

4.53 

40 

1749,0 

0.072 

5.24 

4.93 

3.26 

4,98 

50 

2186.3 

0.067 

3.49 

5.30 

3,51 

5.36 

60 

2623,5 

0.06 J) 

3.70 

5.63 

3.73 

5.69 


^'K = Theraal Conductivity 

D 


TABLE XVII 

THER’vlA.L CONDUCTIVITY OF FIBERS RANDOMLY DISTRIBUTED IN A 3-D SPACE 
FIBER: 0.015" X 0.015" X 0,125” 


ber Voi. 

Vr 



*Kg (IV 

deg-l) 


atio 

Vf 

Q 





(%) 


L ^ 

Al-Ep 

Al-Ep-Ti 

Cu-Ep 

Cu-Ep-Ti 

20 

434.8 

0.115 

2.05 

3.13 

2.06 

3.15 

30 

652 e 2 

0.100 

2.34 

3.57 

2.35 

3.60 

40 

869 0 6 

0.091 

2.57 

3.93 

2,59 

3.95 

50 

1087.0 

0.085 

2.77 

4.22 

2,78 

4.26 

60 

1304,3 

0.080 

2,94 

4.48 

2.96 

4,52 


*Xr, = Thermal Conductivity 




PART V 


SUmkRY OF RESULTS AND CONCLUSIONS 


1. The special technique used in the fabrication of the 
models is very effective for fiber volume ratios of 20 to 50 percent. 
Because of the vacuum technique, very few pits or holes are present in 
the specimens having 20% ratios or lower, and practically none for 
higher ratios. Individual materials can be prepared with reproducible 
properties. 

2. The testing procedure used in this study, i.e,, the "Colora” 
and the ^Comparative Apparatus”, is suitable for evaluation of the thermal 
conductivity of the t\npe of composites used in this study and would seem 
tc indicate that accuracies of better than ^ 3% for the ”Colora” and ^ 2% 
for the ’’Comparative Apparatus” can be obtained. 

3. Effective thermal conductivities were measured for two 
factorial experiments comprising thirty two specimens each. With the 
aluminum fibers embedded in an epoxy matrix, the conductivities ranged 
from 0,95 to 4,61 IV m-1 deg-1. With the copper fibers and epoxy matrix, 
the conductivities ranged from 1.04 to 2,66 K m*^ deg"^. With the aluminum 
fibers and the epoxy-titanium powder matrix, the conductivities ranged from 
1,23 to 12,05 W m“^ deg“^. With the copper fibers and the epoxy-titanium 
powder matrix, the conductivities ranged from 1,79 to 10,05 W m"^ deg"^, 

4. Two statistical analyses; one for the 1/4” specimens, compris 
ing ten samples, and the other, for the 3/4” sample, comprising six samples 
were performed. The fiber volume ratios for the 1/4” samples’ (having the 


epoxy-titanium powder matrix and aluminum fibers) ran^red from 37,5 to 
41.4 percent^ and their measured thermal conductivities ranged from 3,24 
to 4.79 W m“^ dep“^. The fiber volume ratio for the 3/4” samples’ (having 
epoxy-titanium powder matrix and aluminum fibers) ranged from 34 to 35.6 
percent, and their measured thermal conductivities ranged from 5.81 to 
7.02 W m“^ deg"^, For the 3/4” specimens, the mean thermal conductivity 
is 6.25 W m"^ de,e“^, and the standard deviation CT is 0,595. For the 1/4” 
specimens, tlie mean thermal conductivity is 4,27 W m“^ deg"^, and the 
standard deviation CT is 0,534, 

5. Effective thermal conductivities were m;easured for four 
specimens (two 1/4” and two 5/4”) having anisotropic characteristics. 

The specimens v;ere made of copper fibers embedded in an epoxy-’matrix. The 
fiber volume ratios of the four specimens ranged from 39,1 to 42 percent, 
and their measured tliermal conductivities ranged from 9.74 to 15.2 W deg 

6. Effective thermal conductivities were measured for four alumi¬ 
num specimens having the same fiber volume ratios (lOO^’O 3.nd four fiber 
geometries; and four copper specimens having the same fiber volume ratios 
(100%), and four fiber geometries. Their densities ranged from 2,63 to 2,66 
gm cm“^ for the aluminum, and from 8,44 to 8.50 gm cm“^ for copper. Their 
measured conductivities ranged from 8,2 to 17.3 W m"^ deg~^ for aluminum 
specimens and from 10,9 to 13,8 W m*"^ deg“^ for the copper specimens. 

7. The most meaningful factor to be considered for an optimum 
design of the type of composites used in this study are, in the order of 
importance, the conductivity of the matrix, the fiber volume ratio, and the 
fiber cross section as obtained from the factorial experiment. 


137 


8, Anisotropic characteristics for a copper fiber composite 
can be achieved for fiber volume ratios of 40% or hi^^her by just press¬ 
ing the fibers into a predetermined volume. 

9. The thermal conductivity of composite materials, made of 
highly conductive fibers randomly distributed in low conductive matrices, 
is not only a function of the fiber volume ratio, the conductivity of 
the fiber, the conductivity of the matrix as covered by earlier theories, 
but also a function of the fiber geometry and its distribution. Earlier 
theoretical efforts consisted of mathematical models which are highly 
idealized (Maxwell, Rayleigh, Bruggman, etc.); most of these models 
either underestimate or overestimate the thermal conductivities of the 
type of composites presented in this study, 

10, Theoretical equations for the thermal conductivity of a 
matrix filled with fibers, ellipsoidal or cylindrical in shape and random¬ 
ly distributed according to a prescribed fiber volume ratio, have been 
derived. With these equations, a stochastic model has been proposed 
which allows one to simulate realistically composite structures. Limited 
theoretical predictions of the stochastic model seem to check with the 
experimental data, and it shows that additional factors are involved in 
describing composites which are not covered by earlier theories, i.e., 
fiber geometry and its distribution. Theoretical results of com.posite 
thermal conductivity, up to a fiber volume ratio of 15 percent, are 
obtained with the aid of a high speed computer. Higher fiber volume 
ratios can be obtained if a computer with larger memory bank is supplied. 

11. Theoretical predictions of the composite thermal conductivity 
of fiber volume ratios of 20 to 60 percent are obtained by distributing 


fibers in a 3-dimensional space in accordance with a Poisson process, 

This theory was investigated because of the computer limitations of preced¬ 
ing paragraph 10, Experimental results for the fiber geometry of J = 8 

a 

seem to check very well with tliis theory; for ^of 4, 12.4 and 25, this 

a 

theory could be used for fiber volume ratios higher than 40%, 

12, It was found that a specific specimen formulation 
(sample AGED) gives a thermal conductivity similar to that of stainless 
steels but a density only a quarter of that of the steels, 

13. The experimental work indicates the general behavior of 
the thermal conductivity as a function of Kf, Kj^, fiber geometry, and 
its distribution. The data obtained is reliable, and it is recomimended 
to use this data for guidance in developing better, simple, and more 


reliable models. 
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APPENDIX A 

COMPUTER PROGRAM 


The computer program presented in this appendix has been 
designed for fibers with cylindrical cross sections and comprises the 
following tasks: 

1. Generates fibers in space per theory outlined in Part II. 

2. When the desired number of fibers is generated, a "FIBER 
DENSITY CHECK" is made. This is done by dropping a S-line perpendicular 
to the base of the model (Figure 55). The amount of fibrous material 
intersected by this line is calculated as follows: 

a. The distance of a fiber to the -g^-line is first 
determined. This is done by projecting the fiber on to the x-y plane 
and taking the 5^-line as a point P^. Since the coordinates of Pj^ and 
the equation of the fiber in 2-space are known, the perpendicular distance 
d can he easily determined by: 

d = I ^ Byi ^ c I (^.1) 

+ P 

where (xj# coordinates of Pj^ and Ax + By + C is the equation 

of the fiber projection. If d is greater than the radius of the fiber, 
there is no intersection between the Ej-line and the fiber. If d is less 
than the radius, the problem must be examined more carefully. 

b. The coordinates of the intersection must now be found 
and checked to see if they lie within a range slightly larger than the 




fiber projection. If the coordinates lie outside the range, no problem 
exists. If within the range, the following possibilities are to be 


scrutinized. 


(1) If the is contained within the 


calculation of the fibrous material from Figure 35 is: 


E F 

cos 0 


(A-2) 


where EF 


OF = Piber radius. 

OE = Distance of the 2-)-line to the fiber axis. 

P = Angle of the fiber relative to the x-y plane 

Note that the other relations used to obtain EF' are: 


CA-3) 


v;here OD = OF 


and OD* = r 

cos 0 


(A-4) 


(2) If the -Z^-line resides in D or C, a slight variation 


to A may be necessary as can be seen in Figure 35 


r"G = F"E' + E'G 


(A-5) 


using similar triangles, OH D' and EHG: 


E'G _ E'l! 


O’D" O' 


(E’H) X (O'D") 


(A-6) 


where OMl = (fiber radius) X (tan O) 

EMI is a component of the undirected magnitude between the fiber extreme 
and also a component of the 5-line relative to the fiber, 

(3) If the 5^-line is in C, then 

ptTt Qt = pttf Qt . Qf ptl 

The following computer listing is based on the aforementioned outline. 

A detailed discussion on how to use the program is presented in Part IV, 


Section D, 


COMP 1 
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100 COMMOND^DWH^XL^ START»CGUNT 

110 COMMONXM( 500)>XH( 500)>DELTX(500)*C(500)>Z(500)>A(500)>X< 500) 

120 C0MM0NR(500)>S(500)>T(500)>U(500)iSINP(500)>C0SP(500) 

130 DIMENSIONB(500)*Y(500) 

135 C0MM0NDIAG> DIAG2> CNSTF>CNSTM 
140CXL IS THE LENGTH OF THE FIBER 
150CD IS THE DIAMETER OF THE FIBER 
160CD1 IS THE DIAMETER OF THE CYLINDER 
170CH IS THE HEIGHT OF THE CYLINDER 
180C N IS THE SAMPLE SIZE 
190C N1 REP THE NR OF SAMPLES 
200 DATA ST*CT/12345.>0./ 

210 DATA Nl/1/ 

220 DATA XLD.DDD^DID/.125>.015>.750/ 

230 DATA HD/.250/ 

244 CNSTMa.00238 
250 N=0 
260 STARTeST 
270 COUNT=CT 
280 SUMB=0. 

290 SUM1=0. 

300 NADD=0 
310 SUM2=0. 

320 XL=XLD 
330 D=DDD 
340 D1=D1D 
350 H=HD 

360 DIAG=SQRTCXL*XL+D*D) 

370 DIAG2=DIAG/2. 

380 NNP2=2 
390 NNP1=1 

400 CALL FIBERtAC1)>B(1),C(1)>X(1 )jY<1)>Z(1)^SINP(1)>COSP(1)>D£LTX<1)) 
410 CALLEQUAT(A(1)>B(1),C<1)>X(1)*YC1),Z(1),R(I)^S(1)^T(1)^U(1)> 
420&XM(1),XH(1)) 

430 1000PRINT>”SAMPLE SIZE INCREMENT" 

440 INPUT>NN 
450 N=N+NN 
460 NT=0 

470 00201=NNP2»N 

480 2CALLFIBER(A(I),B(I)>C(I)>X<I)#Y(I)^Z(I)* SINPCI)^ COSP<I)> 
490&DELTXCI)) 

500 CALLEQUATCAd),B(I),C(I)>X(I)>Y(I)>Z(I)>R(I)*S<I)»T<I)*U<I)> 
510&XM(I),XHCI)) 

520 J=I-1 

530C P2>P3>P/i ARE COEFF OF P1P2 
540 4 P2=A(I)-A(J) 

550 P3=B(I)-B(J) 

560 P4=C(I)-CCJ) 

570C A1,A2>A3 ARE COEFF OF \J2 
580 A1=A<I)-XCI ) 

590 A2=B(I)-Y(I) 

600 A3 = CCI)-Z( I ) 

610C B1,B2>B3 ARE COEFF OF VI 
620 B1=A(J)-XCJ) 

630 B2=B(J)-YCJ) 

640 B3=C<J)-Z<J) 

650C XN2»XN3#XN4 ARE COEFF OF VIXV2 


COMP 1 


CONTINUED 


660 XNg::^A2*B3-A3*B2 
670 XN3sA3=?^Bi™A5*B3 
680 Xi\'A=Al*B2"A2*Bl 

690C D3 IS THE ABS VAL OF MAG V1XV2 
700 D2"XN2=?^ XN2+ XN 3*XN 3+ XN A*XN 4 
710 D3=SQRT(D2) 

7R0C PERP REP PERP. DIST BETWEEN CENTER LINE OF FIBERS 
730 FEP'P~( P9*XN2+P3*XN3 + P4*XN4) /D3 

7 A n A P £ R P s A B S C p E R P) 

750 IFCAPERP-D)J7^31»31 
760 17 Xf>12sB2*XNA-B3*XN3 
770 XM3~B3«'XN2"B!*XN4 
780 X M 4 -- B ?>' X N 3 •» B 2 X N 2 

790C THE ABOVE REP THE COEFF OF NORMAL TO PLANE V1V3 
800 C isXM2-J^A(J) +XM3*B(J) + XM4*C(J) 

8IOC XM2..XM3^XM4 Cl REP THE COEFF OF PLANE VIV3 

8P0 C2^;A2*ACI )-A!*BCI) 

BOO C3"A3*ACI)-A1*C(I) 

BAOC 01^02^03 REP CONSTANTS OF THE PLANE AND THE LINE INTERSECTING THE 
B5DC PLANE* THE FOLLOWING IS THE SOLUTION OF THE 3 EQ. 

860 D 9 ' A 3 A i * X M 4 + A1 ♦ A1 * X M 2 + A 1 * A 2 * X M 3 

87Q SisCXM3*A1*C2 +A1 *A1*C1 + A1 *C3*XM4)/D9 

BGO S2=<’XM47<A2*C3-A3*C2) -A1*<XM2*C2-A2*C1 ) )/D9 

890 S3~<A3* CXM3*C2 + A1*C1)-C3*(A1*XM2+A2*XM3))/D9 

9C0 KN5=(XN2*PERP)/D3 

9 1 0 XN 6= ( XN 3* PERP ) / D3 

9S0 XN7 = f XN 4* PERP >/D3 

930 S4--'SH-XN5 

940 S5sSg-J-XN6 

950 S6-S3-XN7 

960 IF((ACJ)-S1)*(X(J)-S1))19»19*31 
97C 19 IF(CA(I)-S4)*(X(I)-S4))5*5*31 

9C0 5 NT-NT+! 

990 GOTO2 

1000 31 iFCJ-i)20*20>3 

IGOO 3 J-J“i 

IGSO GO TO 4 

1030 20 CONTINUE 

1040 PRINT*"FOR SAMPLE SIZE =",N," DISCARDED FIBERS’SNT 
1G50 900PRINT*”M0RE FIBERS" 

1060 INPUT*FJ.B 

1070 IFCFIB-!.)1380*790*76 

1080 790PR1NT*"FIBER DENSITY CHECK" 

1090 INPUT*DEN 

HOC I FCDEN-1 . >800*900*900 

1102 800PRINT*"FIBER CONDUCTIVITY CONSTANT" 

IS03 INPUT*CNSTF 

m o PRINT* "VALUES FOR DY AND DX" 

S'jRO INPUT* DY*DX 
1132 NTOT-0 
1134 NT0TL=0 
1136 SUMB"0- 

use syMi" 0 « 

i I 4 0 X C H AN 2 =: D X / 2 . 

USD YCHAM2-;DY/2. 

S I GO YZER0 = 0. 

1 1 70 YZERa=YZ.ER0 + YCHAN2- • 375 
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COMPl CONTINUED 

1190 1 190XZER0 = XCHAN2 

1200 XBOUN=SQRT(•1A0625-YZER0*YZERa) 

1220 1175CALL CALCN(N^XZERO#YZERO.SUMB>SUM 1>NTOT»NTOTL> 

1355 1600XZER0=XZER0+DX 

1360 IF(ABSCXZERO)-XBOUN)1175»1175»1601 
1370 1601XZER0=XCHAN2 
1380 XZERO=XZERO-DX 

1390 1605CALL CALCNCN>XZERG*YZERO*SUMB>SUM 1>NTOT>NTOTL) 

1395 XZER0=XZER0-DX 

1^00 IFCABS(XZERO)-XBOUN)1605»1605*1603 
1407 1603YZERa=YZER0+DY 
1410 IF(YZER0-.375)1190#910#910 
l‘^20 910XNT0T=NT0T 

1425 PRINTl-i26>XNT0T 

1426 1426F0RMAT<IX.12HFIBER HITS =.1P1E12.4) 

1430 XTOTL=NTOTL 

1435 PRINTl^36.XT0TL 

1436 1436F0RMAT(IX.13HMATRIX HITS =.1P1E12.4) 

1437 PRINTl438.SUMB.SUMl 

1438 1438F0RMAT<1X.7HSUMKB =.1F1E12.4.2X.7HSUMKM =.1P1E12.4) 
M40 SUMB=(SUMB+SUM1)/<XNT0T+XT0TL> 

1450 PRINT1885.SUMB 

1460 1885FGRMAT<IX.12HAVERAGE KB =.1P1E12.4) 

1470 SUM2 = SUM2-*-SUMB 

1480 G0TO7 9 0 

1490 1380NNP1=NNP1+NN 

1500 NNP2=NNP1 

1510 GOTO 1000 

1520 76ST0P 

1530 END 

1630 SUBROUTINE FI BERCA.B.C.X.Y.Z.SINP.COSP.DELTX) 

1640 COMMON D.Dl.H.XL.START.COUNT 
1650 CALL RANDCSTART.COUNT.RND) 

1660 2 R=SQRT(RND) 

1670 R=D1/2.*R 

1680 CALL RANDCSTART.COUNT.RND) 

1690 THETA=6.28318*RND 

1700 CALL RANDCSTART.COUNT.RND) 

1710 C=H*RND 

1720 B=R*SIN(THETA) 

1730 A=R*COSCTHETA) 

1750 CALL RANDCSTART.COUNT.RND) 

1760 ALPHA=6.28318*RND 

1770 CALL RANDCSTART.COUNT.RND) 

1780 SINP=RND 

1785 COSP=SQRTC1.-SINP*SINP) 

1790 X=A+XL*COSCALPHA)*COSP 
1800 IFCABSCX)-Dl/2.)7.7.2 
1810 7 Y=B+XL*SINCALPHA)*COSP 
1820 IFCABSCY)-D1/2.)8.8.2 
1840 8Z=C+XL+SINP 
1850 IFCZ-H)9.9.2 
1860 9 IFCZ)2.11.11 

1870 11DELTX=ABSC CD/2.)*SINP*COSCALPHA)) 

1880 IFCDELTX-XL)3.3.4 
1890 4DELTX=XL 
1900 3 RETURN 


COMP 1 


CONTINUED 


19iO END 

1920 SUBROUTINE EQUAT(A>C>X>Y^Z>R>S>T»U>XM#XH) 

i930C THIS ROUTINE YIELDS THE EQ OF PROJ ONTO XY-PLANE IN THE FORM 
!940C Y=MX+H 

1950C THE EQUATION OF CENTER LINE OF FIBER IN THE FORM 
59 60C RX+SY+TZ+U=0. 

1970 XM.= (B-Y)/CA-X) 

1980 XHsB-XM*A 
1990 Ai=X-A 
2000 B1=Y-B 
2010 C!=Z-C 
2020 R=B1+C1 
2030 S=C1-A1 
2040 T=-A1-B1 
2050 C2=A1*B-B1*A 
2060 C3=B1*C-CI*B 
2070 C4=^A!*C"C1*A 
2080 U=C2+C3+C4 
2090 RETURN 
2100 END 

21iO SUBROUTINE RAND(STARTsCOUNT>RND) 

2120 START=10000000./START 

2130C CUT OFF FIRST THREE SIGNIFICANT DIGITS OF QUOTIENT 
21.40C TO GENERATE RANDOM NUMBER. 

2150 IS--START 

2160 FS=IS 

2170 START=START-FS 

21 SO RND=START 

2190 START=START*100000. 

2200 24 COUNT=COUNT+.1 
2210 IFCCOUNT-100000.)4^3>3 
2220 3 COUNT=.l 

22.30C ADD COUNT TO PREVENT REPETITIONS 

2240 4 START=START+COUNT 

2250C CHECK SIZE OF START 

2060 IFCSTART-100000.)5>4^13 

2270 5 1F(START-10000.)6>4>14 

2280 6 IFCSTART-1000.)7»4,12 

2290 7 IF(START-100.)8>4*11 

2300 8 IF<START-10.)9>4*10 

23100 MULTIPLY START BY APPROPRIATE CONSTANT TO OBTAIN NEW START 
2320C WITH FIVE DIGITS BEFORE THE DECIMAL POINT. 

2330 9 START=START*10020.01 
2340 RETURN 

*2350 10 START = START* 1000.2001 
2360 RETURN 

2370 11 START=START*100.002 
23S0 RETURN 

2390 12 START=START*10.00002 
2400 RETURN 

2410 13 START=START/10. 

2420 14 RETURN 
2430 END 

2/-40 SUBROUTINE CAL ON ( N , XZ ERO * YZERO > SUMB> SUM 1 > NTOT » NTOTL) 

2445 COMMON D> D1 H, XL» START* COUNT 

2450 COMMONXMC 500)*XH(500)* DELTX< 500> *C(500)»Z(500> * A< 500)* XC 500) 
2460 C0MM0NRC500)> S(500)* TC 500)*U(500)* SINPC 500)* COSP(500) 


COMPl 


CONTINUED 
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2462 COMMONDIAt>DIAG2>CNSTF»CNSTM 
2465 S'JM=0. 

2490 D040I=1>N 

2500 DIST=ABS(XM(I)*XZERO-YZERO+XH(I))/SQRT<XM(I)*XMCI)+l.) 

2510 IFvDIST-D/2.)41>41>40 
2520 41 SLaPE=XM(I) 

2530 YINT=XH(I) 

2540 XHPRM=YZERO+XZERO/SLOPE 

2550 XVAL2=(XHPRM-YINT)/C SLOPE+1./SLOPE) 

2560 XDELT = DELTX( I) 

2570 IFCCCI)-ZCI))43>43>42 
2580 42 XVAL=X(I)+XDELT 
2590 AVAL=A(I)-XDELT 
2600 GO TO 1182 
2610 43 XVAL=X(I)-XDELT 
2620 AVAL=A(I)+XDELT 

2630 1182 IF((XVAL2-XVAL)*<XVAL2-AVAL>)61»61>40 

2640 61 YVAL=XHPRM-XVAL2/SLOPE 

2650 ZVAL=-(RCI)*XVAL2+S(I)*YVAL+U(1))/T<I) 

2660 DELPX=ABS(ZCI)-ZVAL) 

2670 DELPX=DELPX/SINP(I> 

2680 IF(XDELT-DELPX)72»72*71 
2690 72 DELPX=ABS(C(I>-ZVAL) 

2700 DELPX=ABS(DELPX/SINP(I)) 

2710 71 XLOC=(2./XDELT)*ABS(XZERO-(AVAL+XVAL)/2.)-<XVAL-AVAL)/XDtLT+1. 

2720 MLOC=XLOC 

2730 CA=SQRTCD*D/4.-DIST*DIST) 

2740 CAPRM=ABS(CA/COSP(I)) 

2750 VERT=DELPX*CAPRM/XDELT 

2760 IF(VERT-DIAG2)83.83>84 

2770 84 VERT=DIAG2 

2780 83 IF(ML0C)44^45^46 

2790 44 IFCCAPRM-DIAG2)49»49»48 

2800 48 CAPRM=DIAG 

2810 49 CAPRM=2.*CAPRM 

2820 GO TO 54 

2830 45 IFCVERT-CAPRM)51*51>52 
2840 51 CAPRM=2.*VERT 
' 50 (0T054 

2860 52 CAPRM=CAPRM+VERT 
287 0 GO TO 54 

2880 46 IFCVERT-CAPRM)56.57,55 
2890 56 CAPRM=CAPRM-VERT 
2900 G O TO 5 4 
2910 57 CAPRM=VERT 
2920 G)T054 

2930 55 CAPRM=VERT-CAPRM 
2940 54 SUM = StJK + CAPRM 
2950 40 CONTINUE 
2960 IF(SUM-H)1720^1720#1730 
2980 1730SUM=H 

2990 1720IF(SUM)1670#1670#1660 

3000 1 660S!JMM = SUM/H 

3010 CNSTB=CNSTF*CNSTM/(SUMM*CNSTM+CNSTF*( 1 .-SUMM)) 

3030 200FORMAT<4E12.3) 

3040 NTaT=NTOT+l 
3050 SUMB=SUNB+CNSTB 


COMP I 


CONTI NUEID 


3060 H £ T U RN 

307 0 167 0SUM138 E-3 
3080 NTGTL = NTOTL-^ 1 
3090 RETURN 
4000 END 
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APPENDIX B 

DEFINITIONS OF STATISTICAL AND CHEMICAL TERMS 


1. Population: Any finite or infinite collection of individual things, 
objects, or events. 

2. Sample: A portion of the population. 

3. Distribution: Refers to the distribution of a characteristics or 
a group of characteristics in a population (it is usually obtained 
by measurements). 

4. Randcan Sampling (Simple): This type of sampling is defined by the 

requirement that each individual in the population has an equal chance 
of being the first member of the sample; after the first member is 
selected, each of the remaining individuals in the population has an 

equal chance of being the second member of the sample; and so forth. 

5. Selection of a Random Sample: It is advisable to use a table of 

random numbers as an aid in selecting the sample. 

6. Properties of Distributions: 

Normal or Gaussian distribution, symmetrical bell-shaped 
curve, is completely determined by: 

m, the arithmetic mean (center of gravity) of the 
distribution and 

cr , the standard deviation (radius of gyration about m). 

cr , is the variance of the distribution (second moment about m). 

Z, is the distance from the population mean in units of cr 

and is where x represents any value in the population. 
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P, the proportion of elements in the population which 
have values of Z smaller than any given Z (P is found from table on 
page 549 of reference 28) if Z is known. 

The following graph shows these percentages of the population 
in various intervals of Z. Example: 34.13% of the population will 
have values of Z between 0 and 1 (or 0 and -1). 




7. Estimation of m and (T : 


m = X 


I 

n 



n 

I 

i = 

T 

(n - I ) CTx^ = S (x^) - [ S (x)] 

m and O" will differ from sample to sample. 

8. Confidence Interval: An interval within which we estimate a given 

population parameter to lie (e.g., the population mean ra with respect 
to some characteristic). Inasmuch as estimates of m and CT vary from 
sample to sample, interval estimates of m and 0~ may sometime be 
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preferred to "single-value" estimates. Provided we have a random 
sample from a normal population, we can make interval estimates of 
m or O" with a chosen degree of confidence. 

y - The confidence coefficient is simply the proportion of samples 
of size n for which intervals computed by the prescribed method may 
be expected to bracket m (or O" ). Such intervals are known as confi¬ 
dence intervals. As we would expect, large (size) samples tend to 
give narrower confidence intervals for the same level of confidence. 

9. Statistical Tolerance Limits: For a given population, are limits 
within which we expect a stated proportion of the population to lie 
with respect to some measurable characteristics. 

10, Engineering Tolerance Limits: Are specified outer limits of accepta¬ 
bility with respect to some characteristic usually prescribed by a 
design engineer. 

11, Choice of Null and Alternative Hypotheses: Null hypothesis is con¬ 
sidered to be the hypothesis under test as against a class of alter¬ 
native hypotheses. The null hypothesis acts as a kind of "origin" or 
"base" (in the sense of "base line") from which the alternative 
hypotheses deviate in one way or another to greater and lesser degrees. 

12, Two Kind of Errors: 

Error of the First Kind: If we reject the null hypothesis when it is 
true; e.g., announce a difference which really does not exist. 

Error of the Second Kind: If we fail to reject a null hypothesis when 
it is false. Although we do not know in a given instance whether we 
have made an error of either kind, we can know the probability of 
making either type of error. 


13. Significance Level and Operating Characteristics (OC) Curve of a 
Statistical Test: 

Q. ~ Level of Significance, is the risk of making an error of the 
first kind. ^ - The risk of making an error of the second kind 
varies, as one would expect, with the magnitude of the real differ¬ 
ence and is summarized by the Operating Characteristics (OC) curve 
of the test. Also, the risk increases as the risk ct decreased 
•'Only with large samples can we have our cake and eat it too" - and 
then there is the cost of the test to worry about. 

14. Choice of the Significance Level: The significance level Q of a 
statistical test should be chosen in the light of the attending 
circumstances, including cost. We are occasionally limited in the choice 
of significance level by the availability of necessary tables for 

some statistical tests. Two values of Q , Of = .05 and Ct = .01, 

have been most frequently used in research and development work and are 
given in tabulations of test statistics. We have adopted these "standards 
for the purpose of this study. 

15. Confidence Intervals for the Population Mean When Knowledge of the 
Variability Cannot be Assumed: 

Problem: What is a two sided 100 (1 - a )% confidence interval for 
the true mean m? 

Procedure: 

1. Choose OL . 

2. Compute X and S. 

3. Look up t = tj^ _ J 2/2 n - 1 degrees of freedom (page 

583, reference 28). 
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4. Compute: 

^Upper ~ ^ ^ S/y/n 

\o«er • * - ' 

Conclusion: i.e., we may assert that with 100 (1 - a confidence that 
Xl < m < Xy or the interval from X^^ to Xy is a 95% confidence 
interval for the lot mean; i.e., we may assert with 95% confidence 
that Xy < lot mean < Xy, 

16, Number of Measurements Required to Establish the Mean With Prescribed 
Accuracy: 

Procedure: 

1. Choose d, the alloivable margin of error, and CL , the risk 
that our estimate of n will be off by d or more. 

2. Look up tj^ ^/2 ^ degrees of freedom (page 583, refer¬ 

ence 28). 

3. Compute n or sample size: 



Conclusion: We may conclude that if we now compute the mean X of a 
random sample of size n from the population, we may have 100 (1 - a )% 
confidence that the interval (X - d) to (X + d) will include the 
population mean m. 

17. Amines: Organic derivatives of ammonia; these may be primary, secondary, 
or tertiary, depending on the number of hydrogen atoms that have been 
replaced by organic radicals. 

18. Radical: A group of atoms which occurs in the molecules of a number 
of compoxinds, and which remains unchanged through many chemical 
reactions. Example: Ethyl C 2 HJ (valence 1) 


19. Valence: An intiger representing the nmnber of hydrogen or chlorine 
atoms which one atom of an element can hold in combination. Note: 
groups of atoms, i.e., radicals, also have valences. 

20. Amino Acids: A compound that contains both an amine group and a 
carboxy group. 

21. Carboxyl (or the compound name Carboxylic): A term for the -COOH- 
group, the radical characteristic of all organic acids. 

22. Polyamide Resins: Resins made from polymerized vegetable oil acids. 

23. Polymerized Oils: Vegetable and animal oils which have been heated 
and agitated by a current of air or oxygen. They are partially 
oxidized, deodorized, and polymerized by the treatment and are 
increased in density, viscosity, and drying power. 

24. Polymers: A molecule made up of hundreds, thousands, and even tens 
of thousands of repeating units called monomers. 

25. Monomers: Are usually simple, reactive organic molecules which react 
with one another to link together over and over again and form giant 
molecules, the polymers, which often have molecular weights running 
into the millions. 

26. Polymerization: A union of monomers. 

27. Condensation Products: Some polymers are formed by condensation reac¬ 
tions which, in most instances, involve the elimination of a simple 
molecule such as water or ammonia. 

28. Oxidation: In a broad sense, oxidation is the increase in positive 
valence or decrease in negative valence of any element in a substance. 
On the basis of the electron theory, oxidation is a process in which 
an element loses electrons. In a narrow sense, oxidation means the 
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chemical addition of oxygen to a substance. 

29. Exothermic Reactions: Oiemicals combining to form stable compounds 
and which give off energy in the process. 

30. Exothermic: Refers to a process which is accanpanied by evolution of 
heat. 

31. Epoxy: A prefix in organic nomenclature denoting an oxygen atom joined 
to each of two atoms which are already united in some other way, as 
-C-O-C. 

I_I 

32. "Epon" Resins: Trademark for a series of condensation products of 
epichlorohydrin and bisphenol - A, having excellent adhesion, strength, 
chemical resistance and electrical properties when formulated into 


protective coatings, adhesives, and structural plastics. 

33. Epoxy Resins: Thermosetting resins based on the reactivity of the 

0 

epoxide group _ _ 

34. Epoxidation: The reaction of oxygen with an olifin (a class of 
unsaturated, aliphatic hydrocarbons of the teneral form Cj^ ^ 2 ^ and 
named after the corresponding paraffins by adding "ene" or ylene" to 
the stene: as ethylene, propylene, butenes, etc.) to form an epoxy 


compound. 

35. F (x) = distribution function of the random variable X defined for 

X 

any real number x. 

36. fyCx) = probability density function: f^C^) = ^ 

A 

37. E[X] = the expectation, or mean of X, when it exists is defined by: 

f oo 

E[X] = X f^ (x) dx 


-oc 


APPENDIX C 

YATES' METHOD FOR OBTAINING ESTIMATES OF MAIN EFFECTS AND 
INTEPvACTION FOR TWO-LEVEL FACTORIALS 

YAT ES method FOR OBTAINING ESTIMA TES OF MAIN EFFECTS AND 

_01^ fo r TgQ^LEviEL FAC TOR | ALS 

DIMENSION TABi32o5) «T£mP(32S _ 

_REA D {REA D ER e 500 ) N ? ALPH, T_A_ 

500 for mat ( 15 9 5X 9 ?.F 10» 0) 

_ 

__ 

RNfiNR 

read{READER*505)((TA8(I♦J) 9 J«1»3> *I# 1 ,nR) 


505 

FORMAT{A5« 5X« 

2F10.0) 

C 

PUT RESPONSE 

IN SAVE area 

DO i 1^1.NR 

1 

TEMP< I )gfTAB< I 

.3) 

C 

PERFORM 

ADDITIONS AND SUBTRACTIONS TO OBTAIN GT.Ga.ETCo 

c 

FINAL VALUES 

OF G ARE IN COLUMN JJ OF TAB 

MRH?^NR/2 


.03 





DO 4 KgiJ.N 


IE#2 



lOFI 



DO 2 IFlsNRH 

TAgC IjjJJS«TAB{ IE* J>J^TAB< 10. J) 

TAS<I+NRH.JJ)«TAB<IE.J)-TA8(lOeJ) 


2 jOai 
IDI 


J«JJ 
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4 JjaiDUM _^_ 

__*J 

C divide values of G by 2»»(N"1) _ 

DO 6 I« ltNR___ 

6 TAB<I*5)tfTAB<ItJJ)/(2«»»(N»l)) _ 

c calculate S2 which is the sum of the squares Of g s _ 

_ S240. ____ ___ 

DO 10 I#ltNR _ 

_ IF(TAB(I»2)~3«)10t8»8 _ 

6 s2 FTAB(I>JJ)*TAB(ItJJ)»S2 _ 

10 CONTINUE __ 

c calculate nu and W _ 

_ UN4RN~0»5»(AN«AIM-»AN’*-2« ) _^_ 

S24S2/(RN*UN) ______ 

W<^TA*SQRT(RN*S2) __ 

_ pO 11 ItfltNR _ 

11 wR I TE(PRINTERt510)TAa(ltl)tTAB(I«2> «TEMP(I)«TAb(1<JJ)*TAB(I i5) 

510 rORMAT( 10X.A5t5XtEl4.2«3E20«8> __ 

_ F£ED(3) __ 

_ W^lT E<PRINTERt515)S2»W __ 

515 format (IPX ♦3HS2«F15»at IPX t2HW(WFlQ»5»//) _ 

C_ IF ABS(G(I>) IS greater THAN M THEN PRINT G(l) _ 

_ DO 14 1«1«NR __ 

_ IF(AB5(TAB(I*JJ))»W)14,14tl2 __ 

_1 2 WRlTE(PRlNTERt520)TA8( I ♦ JJ) ___ 

520 F 0RMA T<E20.8) _ 

14 continue _ 

end _ 
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APPENDIX D 


DISTRIBUTION OF THE NEAREST NEIGHBOR IN A POISSON DISTRIBUTION 

OF FIBERS IN SPACE 


Reference 40 was used for the following analysis. If particles 
or fibers are distributed in a S-dimensional space in accordance with a 
Poisson process at a mean rate or intensity 1/ per unit volume, then 


p (K; y V) =e"^^ 

K! 


(D-1) 


considers the probability of finding exactly K particles in a fixed 
volume V, In particular, the probability of no particles in a volume V 
is 


p (0. yv) = e 


-yv 


(D-2) 


and the probability of one or more particles is therefore 


1 - e 


- yv 


(D-3) 


From Appendix B, definition 35, 


(X) = P X<,x*| 


and from definition 36 


(D-4) 


fy (X) = 


dx 


Px Cx) 


(D-5) 


If S is the distance between an individual fiber and its nearest 
neighbor, then 
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4 „ r 




( 0 - 6 ) 


f^(x) 


d_ ^ 


4T7i/X" exp U i-ni/x" 


(0-7) 


To find the mean distance 


. I’ [s]. 


between a fiber and its nearest neighbor. 


definition 37 of Appendix 0 is used; 


[s] = 


\ X f,., (X) dx 


\ ■' ^(« 


( 0 - 8 ) 


hence 


[s] = 


J_ ■) r 

4 77 I I 


j = 0.554 1 / ‘ V 3 


( 0 - 0 ) 


:here is equal to the number of fibers per unit volume. The fiber 


volume ratio 




(D-IO) 


Substitutinr eauation (D-10) into (D-9) for , the following expression, 


for S , is obtained: 
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where Vj^ is the fiber volume ratio and is the volume of 1 fiber. Once 
8 is determined, the average thermal conductivity is found in a similar 
procedure outlined in Part II, Section A. 

Or 


Kb = 


Kf 

§)k, 


(D-12) 


where L is the length of the model. Note that if 8 is equal to zero, K 


is equal to K^. Results from this theory are presented in Tables XIV 


through XVII. 
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